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CHAPTER I 


INTRODUCTION 


1.1 General 

A new and comprehensive view of Earth is needed to cope with 
environmental problems, as well as with difficulties caused by an 
expanding population and the depletion of natural resources. Earth 
must be viewed in its entirety as problems with the air-ocean-land 
system are global in extent. Earth-orbiting satellites, including the 
Earth Resources Technology Satellite (ERTS), have been developed which 
are able to remotely sense the Earth features as well as the data 
collected by specifically designed data collection platforms (DCP).^ 

^•2 The Alabama ERTS Project 

The University of Alabama, with the participation of the Geological 
Survey of Alabama and the Marshall Space Flight Center, has undertaken 
a study of the feasibility of applying remotely sensed data to the 
management of natural resources and to the improvement of environmental 
quality in Alabama. The accomplishment, of this purpose will require the 
following: 

(1) Identification and education of users with regard to the 
potential benefits to be derived from space-acquired data, 

(2) Timely interpretation and dissemination of beneficial infor- 
mation to ultimate users, especially policy makers and regula- 
tory agencies. 
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(3) Analysis of whether or not information from remotely sensed 
data results in a significant improvement of the user's 
decision-making ability and actions related to management 
of natural resources and environmental quality . ^ 

The overall objectives of the Alabama ERTS project are summarized 
in, outline form in Table I. In addition. Tables II and III show, respec- 
tively, areas of application of ERTS data and professional people who 
could derive benefits from ERTS data. 

Because water is one of the major natural resources and plays an 
important role in the development of other resources, particular empha- 
sis is given to water resources both as to quantity and quality in the 
Alabama. ERTS project. Therefore, one of the primary purposes of the 
environmental phase of this project is to test the feasibility of using 
remotely sensed data in conjunction with ground truth data to monitor, 
predict, and manage water quantity and quality in our waterways. 

There are several applications for this concept to the monitoring 
and management of environmental quality. One application is concerned 
with the fast and accurate acquisition of data necessary to meet the 
needs of water resource managers, which include both regulatory agencies 
and private industries.^ Results from data collection platforms stra- 
tegically placed in a water basin, together with satellite imagery, may 
be used to characterize environmental factors and provide means to record 
changes in environmental quality. In addition, the ERTS data from the 
DCP may be employed in monitoring environmental changes with regard to 
enforcing state and federal regulations. The collected data may be 

stored in an information system which would allow access to specific 

2 

data according to . the needs of the user. 



TABLE I 


(BJECnVES OF ALABAf'R ERTS PROJECT^ 

A. TO DETEmi'C THE APPLICABILITY OF RBCTELY S0KED DATA 
FROn ERTS FOR; 

(1) LAND USE 

(2) IN'filTORY AND fRN.A®ENT OF NATURAL RESOURCES 

(3) IffROTeENT OF THE QUALITY OF BWIRaifNT 

B. TO DISSB1INATE INFORmTiai IN fURfB IttST SUITfflLE FOR 
ULTimTE USERS: 

(1) PUBLIC POLIQ' TECIWCIANS 

(2) DECISIffl IRKERS 

(3) PRIVATE INDIBTRIES 


(A) PRIWTE CITI2BIS 



TABl£ II 


USES OF ERTS mTA^ 


FLOOD CaiTROL 

DISASTER DETECTION 

SOIL STUDIES 

DAmCE EVALUATiai 

RESOURCE INVENTORY 

SEDIf'ENT TRANSPORT 

SURFACE HATER SUIDIES 

TRAFIC STUDIES 

HlfERAL EYPLDRATIffI 

EROSIdl CONTROL 

GROUND HATER STUDIES 

IRRIGATION 

1 

CD 

WATER TePERATURE STUDIES 

GRam TRGIDS 

cmp coTDiTiais 

RECEEA.TICN 

SURVEYING m mFING 

WNAGBeiT 

AIR QUALITY' mummi 

PESTICIE STT5T1ES 

WATER QUAUPi' f-TOGEIBTT 


^^S^N fftt) REGiaiAL PLANNING 



TABLE III 


POTENTIAL USEfS OF ERTS DATA^ 


mm PLANNERS 

CIVIL ENGirCERS 

REGia^L PLANNERS 

CH01ICAL ENGINERS 

FORESTERS 

AGRICULTURAL B'lGINERS 

EOIDGISTS 

MINING ENGIItERS 

ECOU3GISTS 

GEOGRARIERS 

HYDROLOGISTS 

UPNOU3GISTS 

AGRO^Of'IISTS 

ENTOMOLOGISTS 

BIOLDGISTS 

ARCHITECIS 

PHYSICISTS 

ARCHAEOLOGISTS 

ASTRaWERS 

DB«RAPHERS 

aeiisTS 

WTYERS 

AGRICULTURISTS 

UNIVERSITY FACULTY IBBERS 
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In order to monitor a water basin in an effective and optimum 
manner, it is necessary that the locations of the data collection plat- 
forms be selected such that measurements of water quality parameters at 
these locations will be indicative of water quality over the entire 
basin. This study is concerned with developing methodology for selec- 
tion of these optimum monitoring sites as specifically noted in the 
ob j ec tives given below , 


1.3 Ob i actives 

The objectives of this study were: 

(1) to develop the methodology, based upon a systematic investi- 
gation, for selection of desirable locations for placing 
remote sensing devices (data collection platforms) in a 
waterway such that the data collected will be indicative 

of the water quality over the entire basin. 

(2) to select desirable locations for placing data collection 
platforms in a portion of the Black Warrior River basin. 

(3) to develop a mathematical model which may be employed to 
predict water quality in the Black Warrior basin. 



CHAPTER II 


MONITORING SYSTEMS 

2*1 Conventional Automated 
Monitors 

Changes in water quality, such as those caused by storms, industrial 
spills, and flow changes from impoundments, often occur suddenly and 
affect the concentration of many substances of particular interest. For 
this reason, continuous monitoring of water quality is advantageous in 
situations where these abrupt changes are likely to occur. Because of 
the manpower and time requirement involved in standard manual sampling 
and analysis, testing in water pollution control work has been moving 
more toward instrumentation. 

A typical installation utilizing a continuous, automated monitor is 
shown in Figure 1. Most automatic monitors have three basic components — 
sensor system, analyzer phase, and output phase. The sensing element, 
which is the part of the system in contact with the sample, may be either 
immersed in the stream or set in flow cells through which the sample is 
pumped. Because of advantages in design and maintenance, most systems 
use the flow cell principle; however, when the sensor is directly in the 
stream, the sample is not affected by pumping, temperature changes, or 
time of travel through the instrument. 

The types of sensors currently available fall into several cate- 
gories; (1) electrochemical sensors, (2) sensors based on colorimetric or 
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Figure 1. Conventional Monitoring Inatallation^ 
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light-scattering measurements, (3) sensors measuring temperature through 
use of a thermocouple, and (4) sensors measuring physical parameters 
such as velocity. 

The analyzer phase of the monitoring instrument converts the signal 
from the sensor into a voltage to apply to the output phase. The analy- 
zer may be designed to receive signals from one or more than one sensing 
element. 

The output phase of the instrument presents the measured value in 
the necessary units — pH units, milligrams per liter, micromhos, etc. — 
and records it permanently. This component normally has a meter panel 
on the face of the instrument to indicate output to the recording 
devices . 

Where a number of monitoring devices are interconnected to provide 
a simultaneous evaluation of water quality in a river system, a tele- 
metry system to provide remote handling of the data has almost become a 
requirement due to the need for speedy collection and analysis of the 
data. ^ 

2. 2 Description of ERTS Data 
Collection Platforms 

The instrumentation used in the ERTS program for monitoring water 
quality, xdiile in some respects similar to conventional instrumentation, 
embodies some innovative features in its design and construction. 

Figure 2 shows a sketch of one of the data collection platforms to 
be used in the Alabama ERTS project. Its primary components consist of 
the radio transmitter inside the upper housing, the antenna located on 
top of the housing, and the sensors, which are included in the lower 
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PRIMARY COMPONENTS: 

1 - Metal Pole (support) 

2 - Sliding Fastener 

3 - Housing for Analyzer Unit 

Battery, and Electronic 
Unit 

4 - Antenna 

5 - Connector from Sensor to 

Analyzer Unit in Housing 

6 - Quick Release Pin 

7 - Sensor 



Figure 2. ERTS Data Collection Platform 
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extremity of the DCP. The instrumentation of the DCP will be fastened 
to a metal pole which has been driven securely into the river bottom. 

The fastener will allow vertical movement :of the instrument along the 
pole, facilitating its removal for recalibration and other servicing. 

The DCP will be submerged except for the antenna and a small portion of 
the housing, allowing monitoring at a five-foot depth. Another attrac- 
tive feature of the sliding fastener mentioned above is that it will 
allow the DCP to float and maintain the five-foot monitoring depth 
regardless of changes in the water level. 

The DCP will be one of three basic units in the overall data collec- 
tion system, the other two being the ERTS satellite and the ground receiv 
ing station. As now planned, the DCP will be activated by a timer 
shortly before the satellite is due to pass over in any of several adjac- 
ent orbits within a line of sight of the DCP. The DCP will remain opera- 
tional until the satellite has passed overhead, after which the instru^ 
ment will be deactivated. The water quality measurements taken during 
this period of operation are immediately beamed to the satellite by 
means of a radio transmitter and antenna which are part of the DCP 
itself. Computers on board the satellite will correlate this data 
received from the DCPs with imagery taken at the same time by sophisti- 
cated scanners on board the satellite and transmit the data to the receiv 
ing stations on the ground.^ 

Two requirements regarding installation of the DCP are that it be 
(1) near the shore, where possible damage by contact with river traffic 
in the main channel, will be minimized, and (2) in a position such that 
the line of sight to the satellite may be maintained. This latter 
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requirement would exclude areas having high or overhanging cliffs on 
either bank of the river. 

^•3 Comparison of the ERTS Data 
Collection Platforms with 
Conventional Automated 
Monitors 

In describing both the conventional automated monitor and the DCP 
to be used in the ERTS program, several advantages of using the DCP 
system became apparent. The primary advantage of the DCP is that it 
allows immediate correlation of water quality data with imagery taken 
from the ERTS satellite, in order that a methodology might be developed 
for detecting the occurrence of significant changes in the quality of 
water. Another advantage is that this DCP data could be used for develop^ 
ing means of detecting changes in basin characteristics and the constit- 
uents of runoff from the ERTS imagery. 

Another attractive feature of the DCP is that after the satellite 
comes down there is a possiblity that the DCP could be utilized in the 
same manner as the conventional automated monitor. 

2.4 Parameters Monitored 

In the selection of water quality parameters to be monitored, it is 
important to select those which would be of the most benefit to the ulti- 
mate users. It was decided that in the Warrior River basin the parame- 
ters dissolved oxygen, temperature, pH, and specific conductance would 
provide the most beneficial information for the management of basin water 

o 

quality. Table IV, which lists parameters according to frequency of 
usage in state water quality standards, indicates that, with the excep- 
tion of conductance, which is used less than 20 percent of the time, the 
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TABLE IV. 

FREQUENCY OF PARAMETER USAGE IN WATER QUALITY CRITERIA 
OF STATE STANDARDS 6 


Uniform 

Frequent 

Infrequent 

Rare 

(100%) 

(99-50%) 

(49-20%) 

(19-0%) 

DO 

Radioactivity 

Arsenic 

Bottom Deposits 

pH 

Public Health 
Service Drink- 

Barium 

Chromium (+3) 

Coliform 

ing Water Stds. 

Cadmium 

Electrical 




Conductance 

Temperature 

Total Dissolved 
Solids 

Chromium (+6) 

Ammonia 

Floating Solids 


Fluoride 

Acidity 

(Oil-Grease) 


Lead 

Alkalinity 

Settleable Solids 


Selenium 

CCE 

Turbidity and/or 

Color 

Silver 

Hydrogen Sulfide 



Suspended 

Solids 

Pesticides 

Taste-Odor 


Turbidity 

Sodium 

Toxic Substances 


Chloride 

Iron 



Copper 

Plankton 



Nitrate 

Foaming 

Substances 



Phenols 

Boron 



Phosphate 

Manganese 



Sulfate 

Hardness 



Color 

BOD 



Cyanide 

MBAS 


Zinc 
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selected parameters are those uniformly used in formulating state water 
quality criteria. In view of the fact that state regulatory and plan- 
ning agencies are anticipated to be one of the principal users of ERTS 
water quality data, the choice of these four parameters appears to be 
advantageous . 



CHAPTER III 


SELECTION OF THE WATER BASIN AREA 

3.1 Background 

The original plans for this study called for the use of ten DCPs in 
monitoring the water quality of several river basins on a statewide basis 
Later developments indicated that the number of available DCPs would be 
reduced to the five now allotted to this portion of the Alabama ERTS 
investigation. Based on this reduction, it was felt that the investiga- 
tion should be confined to a single river basin to avoid a thin coverage 
of a multiple basin area, where the data collected would be from rela- 
tively isolated and unrelated points. Also, simultaneous evaluation of 
water quality with ERTS imagery over a large geographical area would be 
less significant than that produced for a smaller and more extensively 
monitored area.® 

3-2 River Basin Selected 
for Investigation 

The sector of the Black Warrior River from river mile 385.0 to river 
mile 335.0 was chosen for evaluation in testing the utility of the DC? 
concept for monitoring water quality. Shown in Figure 3, this stretch 
of the river extends from the confluence of the Locust and Mulberry forks 
downstream to a point approximately three miles below Oliver Lock and 
Dam. These limits were selected, in part, because they cover all sig- 
nificant pollutional effluents received by the Black Warrior River from 
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the highly industrialized area in the vicinity of Tuscaloosa. Another 
unique factor concerning this sector is that the relative quality of the 
water both upstream and downstream of these limits is good while portions 
of the area considered are grossly polluted, especially during the low 
flow periods occurring in the summer months . 

The reasons for choosing the Warrior River Basin for the study are 
summarized as follows: 

(1) overall convenience of the basin in regard to data collection 
and maintenance of DCPs to the primary investigating groups 
at The University of Alabama, the Geological Survey of Ala- 
bama, and Marshall Space Flight Center. 

(2) relatively unrestricted use of the river as a receiving body 
for disposal of sewage and industrial wastes. 

(3) location on the river of several major industries, such as 
chemical, petrochemical, iron products, coking, pulp and 
paper, and asphalt operations, which discharge wastewaters 
with varied constitutents . 

(4) location of several dams and hydroelectric power generation 
facilities at Bankhead, Holt, and Oliver, which would allow 
the investigation of their effects on the quality of water in 
an imp oun dmen t . 

(5) presence of both stratified and unstratified conditions in 
the reservoirs, primarily due to the deep, slow-moving waters 
in Bankhead and Holt pools and the shallow, swifter waters in 
Oliver and Warrior pools. 
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(6) presence of both relatively unpolluted waters, such as 
Bankhead Pool, and grossly polluted areas, allowing com- 
parative modeling techniques various degrees of water 
quality, 

(7) potential for additional investigations, such as evaluating 
the moderating effect of a clean stream (Mulberry Fork) join- 
ing a relatively polluted stream (Locust Fork), 

The above reasons for selecting the Black Warrior River point out 
the numerous advantages inherent in selection of this area for investiga- 
tion during this part of the environmental phase of the Alabama ERTS 


project. 



CHAPTER IV 


REVIEW OF LITERATURE 

Since the ERTS remote-sensing concept involving correlations of 
imagery taken from the satellite with water quality data obtained by the 
DCPs is new and unprecedented, reliance must be placed on previously 
reported methods for selection of DCP locations and for the development 
of stream quality models. 

4 . 1 Dissolved Oxygen Models for 

Flowing Streams and Impoundments 

The first dissolved oxygen model for predicting oxygen balance in 

Q 

a flowing stream was developed by Streeter and Phelps in 1925. The 
formulas, which are based upon two velocity constants, describe the oxy- 
gen balance in a stream as a function of distance (or time) from a waste 
load discharge point. These two parameters are kj^, the deoxygenation 
velocity constant, and k 2 , the reaeration velocity constant. These con- 
stants describe the activity in the stream in an all-encompassing fash- 
ion, with the effects of several known interacting factors such as 
photosynthesis and bottom deposits considered as being included in the 
k^ and k2 values. 

Other investigators have proposed models which differ somewhat, but 

all employ the original S treeter-Phelps formulation as the basis of their 

work. Goodman^^ was one of these subsequent investigators who developed 

a mathematical model to apply to flowing streams containing no reservoirs 

and not influenced by estuaries. This particular effort toward water 

19 
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quality modeling was based on modifications of the original Streeter- 
Phelps equation and considered dissolved oxygen concentration as the 
principal criterion of stream quality. These equations determining 
changes in BOD and DO require input values for deoxygenation and reaera- 
tion velocity constants, settling out of BOD to bottom deposits, resus- 
pension of BOD from bottom deposits, and oxygenation by photosynthetic 
processes. 

More recent work involving modeling of flowing streams has been 
reported by the Texas Water Development Board. The result of this work 
was the development of a computer program called QUAL-I^^ that is capable 
of producing a time history and spatial distribution of not only BOD and 
DO but also temperature and as many as three minerals. This is accom- 
plished within the framework of a completely-mixed, branching stream or 

canal system with multiple waste inputs and withdrawals. In addition, 

12 

this agency has developed the DOSAG-I program, which is used to simu- 
late the spatial and temporal variations in BOD and DO under various 
conditions of temperature and headwater flow. 

One of the more useful recent studies done on the subject of model- 
ing dissolved oxygen was conducted by Pyatt, who developed equations 
for predicting organic matter and DO deficit not only for a free-flowing 
stream but also for an impounded waterway. These formulations are based 
on the original Streeter-Phelps equations with added "error” terms to 
account for the usually omitted factors such as bottom deposits and 
photosynthesis which affect the deoxygenation and reaeration rates in a 
stream environment. 

The modeling of the reservoir-type situation is particularly note- 
worthy since relatively few investigators have even attempted such a 
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study. Due to the lack of previous knowledge in this particular area, 

complete mixing was assumed for the reservoirs, although this is not 

actually the case in most impoundments. 

14 

Churchill and Nicholas also investigated the changes occurring in 
the quality of water during its passage through Tennessee River reser- 
voirs and during long storage in impoundments. However, these studies 
were primarily observations in which no attempt was made toward develop- 
ing generalized equations for predicting changes in impoundments. 

^•2 Water Quality Monitoring 
Systems 

Quite a number of investigating groups have undertaken to set up 
systems to monitor and control water pollution. Describing all of these 
in detail would be a formidable task; therefore, only the ones considered 
applicable to this study are reviewed. 

McCormack and Perlis^^ developed a method of optimizing the number 
and locations of measurement stations needed for a particular monitoring 
program. In developing this procedure they applied theories of system 
optimization to a polluted stream model in which the primary dependent 
variable was dissolved oxygen. The stream models were subject to random 
variations and environmental changes. Measurement error was a function 
of the number and position of the measurements, the sample size, and the 
time between measurements. The developed policy minimized an integral- 
type function involving mean square error and actual measurement cost. 

Other studies, which may not be directly concerned with optimizing 
locations of monitors, are strongly related to the use of water quality 
data in water resources management and therefore pertinent to the entire 
ERTS project. A Harvard research teara^^ on the Lehigh River was the 


O 
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first group to conduct regional water resource management studies using 
systems analysis techniques. In recent years, both federal agencies and 
universities have used systems analysis on river basins, including water 
quality studies on both the Columbia and the Delaware river systems. 

Another group to enter this field of basin management was the Ohio 
River Valley Water Sanitation Commission (ORSANCO)^^ which in 1960 built 
an automatic field monitor and a central receiving station to aid in 
management of water quality in the river basin. Later, more field moni- 
tors and data processing facilities \7ere added to make up the integrated 

system now in operation. 

18 

Testerman reported on a system for recording and transmitting 
digital data at a remote water quality monitoring station that is unat- 
tended. Signals from transducers, which measure water quality charac- 
teristics, are converted to digital signals and recorded on magnetic 
tape. This unit can be contacted from a central unit for playback of 
the day’s recording with the transmitted data being recorded by teletype 
at the central station. The techniques developed in this study hold 
promise for making water quality measurements in remote areas and at 

numerous sites and reporting to one central station. 

6 

The United States Army Corps of Engineers has also conducted an 
Interesting study involving the control of impoundments to augment flow 
based on information reported by a monitoring station concerning dilution 
requirements . 



CHAPTER V 


DEVELOPMENT OF THE MODEL 


5, 1 General 

The steps outlined in this chapter were undertaken to develop a 
technique which would ultimately lead to the selection of the desirable 
sites for placing the ERTS data collection platforms. It was intended 
that the methodology developed in this study would find general applica- 
tion for the placement of remotely-sensed monitoring stations^ 

Since dissolved oxygen is normally considered to be the most impor- 
tant parameter in defining water quality, it was decided that the pro- 
cedure for selecting DCP monitoring sites would be oriented toward locat- 
ing the critical dissolved oxygen concentrations in the river basin. 

For this reason, the modeling technique considered in this study was 
centered around simulating the dissolved oxygen concentrations normally 
found in the river. 

While dissolved oxygen is considered to be the principal parameter 
of interest, the additional parameters of pH, temperature, and conduc- 
tivity are also to be monitored. Factors" such as mine drainage and run- 
off from agricultural lands, which would affect these latter three 
parameters, were also considered in developing the methodology for 
selecting the monitoring sites. 
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In deciding to select locations of monitoring sites based on the 
most sensitive points for dissolved oxygen, it was felt that the major 
effects on the other three parameters could also be detected at these 
same locations. The opposite of this, that is, selecting locations 
based primarily on monitoring pH, temperature, and conductivity would 
not necessarily produce locations of meaningful dissolved oxygen con- 
centrations. The reasoning here was that while dissolved oxygen content 
is very much dependent on time (and therefore distance) , factors that 
would be indicative of changes in pH, temperature, and conductivity, 
such as dissolved solids content, are largely independent of time of 
flow, depending instead on the amount of dilution received. 

For the reasons discussed above, the technique developed in this 
study does not simulate pH, temperature, or conductivity, but it does 
include the capability of receiving inputs for factors affecting these 
parameters should it become desirable to include these values in model- 
ing studies. 

Factors Affecting Oxygen 
Balance in Streams 

The factors responsible for the occurrence of oxygen variations in 
streams as well as for the magnitudes of these variations are listed by 
Goodman^® as: (1) deoxygenation, (2) reaeration, (3) algal activity, 

(4) benthal demand, (5) settling and resuspension of BOD, (6) temperature, 
(7) streamflow, and (8) sunlight. The factors which serve to effect a 
decrease in dissolved oxygen concentration are deoxygenation due to 
oxidation of organic matter, benthal demand (bottom deposits) , resuspen- 
sion of organic matter having an oxygen demand, an increase in tempera- 
ture of the water, and a decrease in streamflow. Conversely, those 
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factors which would contribute toward an increase in the dissolved oxygen 
are reaeration, the photosynthetic activity of algae in combination with 
sunlight, settling of organic matter having an oxygen demand, a decrease 
in temperature, and the diluting effect of an increase in streamflow. 

In Figure 4 a diagram is presented depicting activities which would exert 
some measure of influence on the BOD and DO systems of the total river 
environment. 


5 . 3 Streeter-Phelps Equation 

The first attempt to mathematically relate the factors affecting 
the oxygen balance in a natural stream was performed in 1925 by Streeter 
and Phelps.^ These investigators combined all the factors adversely 
affecting dissolved oxygen into a basic deoxygenation equation, while 
those factors contributing to an increase in dissolved oxygen were com- 
bined into a basic reaeration equation. These two equations are shown 
below in both their differential and their integrated forms. 


f - -hi- 

L(t) = (la) 


dt 


-IC2D 


D(t) = D^(e"^2t) 


( 2 ) 

(2a) 


In the above equations L(t) and D(t) are, respectively, the concen- 
tration of organic matter remaining at time t and the dissolved oxygen 
deficit at time t. The terms La and Da are, respectively, the concentra- 
tion of organic matter and the dissolved oxygen deficit at t = 0, while 
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and k .2 are, respectively, the overall deoxygenation and reaeration 
rate constants. 

Equations (la) and (2a) were combined to give the well-known dis- 
solved oxygen sag equation shown below in Equation (3) . 

D(t) = (e-klt + Da(e-^2t) (3) 

k2“ki 

This equation, developed for the case of a free-flowing, non-impounded 
stream, has remained the basis of nearly all water quality research 
related to dissolved oxygen. 

5.4 Modifications of the Original 
Streeter-Phelps Equation by 
Pyatt 

While the original dissolved oxygen sag equation developed by 

Streeter and Phelps has remained largely unchanged, some investigators 
1 q 

such as Pyatt have attempted to make refinements, particularly in the 
area of altering the rate constants to account for factors such as photo 
synthesis and bottom deposits in some manner other than grossly combin- 
ing them all together in the kj^ and k 2 rate constants. 

Pyatt has proposed accounting for the additional factors in the 
form of certain "error’* terms. He has added a factor r, the deoxygena- 
tion error term, to the general deoxygenation equation to obtain: 

ft » -k,L + r (4) 

which, when properly integrated, gives: 


L(t) = (La - 


(5) 
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Similarly, Pyatt has added a factor s,. called the reoxygenation error 
term, to the dissolved oxygen sag equation to obtain: 


dP 

dt 


kjL - k. 2 p + s 


( 6 ) 


which, when integrated gives: 


D(t) 


k2-ki 


r 

k2"ki 




+ 1_ (r+s)(l-e ^2^) + Da 
^2 


(7) 


The error terms r and s are usually set equal to 0,005. All other 

variables in Equations (4) through *(7) are as previously defined. 

One drawback to the use of Equations (5) and (7) is that they are 
limited to use with a flowing stream. Recognizing this limitation, as 
well as the fact that most streams are now regulated to some extent, 
Pyatt sought to develop equations for the concentration of organic mat- 
ter and the DO deficit in an impoundment. His efforts resulted in the 
following equations : 


“At 

L(t) = e At + ^ 


(1-e-At) 


(8) 


D(t) = 


Do - 


ZDin ZL,- „ 

+ -A- do - -^) 


B A-B 


-Bt 


+ 


^(L - 

B-A o A 


e~At + ^ (D 

B 


in 


+ MiH) 


(9) 


where the variables are defined as follows: 


Lq = concentration of organic matter in the reservoir before the 
input of waste flow from upstream 
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= dissolved oxygen deficit in the reservoir before the input 
of waste flow from upstream 

Lin = incoming concentration of organic matter 

Din incoming dissolved oxygen deficit 

W = outflow -r volume 

Z = inflow’ V volume 

IC 3 = deoxygenation rate constant for sludge deposits on bottom 
of reservoir 

A = k-i 4* IC 3 + W 

B = k2 + W 

One limitation of Equations ( 8 ) and (9) is that complete mixing in 
the reservoir must be assumed, although this is often not the actual 
case* Even with this built-in source of error, which would probably 
cause predicted values of dissolved oxygen concentrations to be slightly 
lox^er than those actually found, these equations have been found to be 
particularly useful in simulating the dissolved oxygen profile of a 
stream. Pyatt utilized these equations in simulating the DO profile in 
an actual river basin containing impoundments, showing that the formulae 
developed in his study do provide accurate results.. 

5.5 Application of Formulas 
to Warrior River 

I^Jhen considering the stretch of the Warrior River with which this 
phase of the Alabama ERTS project is concerned, it became apparent that 
no single equation can be applied over the entire run of the river from 
the confluence of the Locust and Mulberry forks to several miles below 
Oliver Lock and Dam. 

This area of the river, which includes all of the Holt and Oliver 
pools and portions of the Bankhead and Warrior pools, is, in the 
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strictest sense, impounded over its entire length, suggesting that 
Equations (8) and (9), developed for use with impoundments, should be 
used for this entire stretch of the river. This would certainly hold 
true for the Bankhead and Holt pools, which, in addition to containing 
the deep, slow-moving waters characteristic of impounded streams, are 
also regulated entirely by the release of water from Bankhead and Holt 
dams. However, in the case of Oliver Pool and the portion of the Warrior 
Pool included in this study, the waters are shallower and faster-moving 
than in Holt and Bankhead pools upstream. Another factor to be con- 
sidered is that Oliver Dam is a free-flow dam with the water flowing 
over the crest. As a result, these waters are subjected to much less 
regulation than that found at Bankhead and: Holt dams, which are not of 
the free-flow type. 

For the reasons enumerated above, it was felt that Bankhead and 
Holt pools are impoundments requiring the use of Equations (8) and (9) 
for dissolved oxygen simulation. However, the Oliver and Warrior pools, 
though technically impoundments, more nearly approach the conditions of 
a flowing stream, and warrant the use of Equations (5) and (7) which 
were developed for free-flowing conditions. 

It should be emphasized that this step is of paramount importance, 
requiring careful consideration of the stream under investigation before 
selecting the appropriate equations for use in simulating the dissolved 
oxygen profile. 

Reach Concept Applied 
to the Model 

Due to an absence of industries and other such inputs along the 
impoundments modeled by Pyatt, each of the impounded areas was considered 
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as one reach or segment having no subdivisions. This approach lacked 
sufficient detail to model a complex basin such as the Warrior; however, 
the modeling equations developed by Pyatt can be applied to the Warrior 
River by subdividing each impoundment into smaller reaches that are 
capable of describing additional inputs. 

This concept essentially involves separating the entire length of 
the river under consideration into a number of smaller reaches, such 
that each reach contains only one input, if any, and also to insure that 
the physical characteristics of the river remain the same throughout 
each individual reach. 

In order to follow the above guidelines, a new reach is begun when- 
ever one of the following events or structures is encountered: 

(1) entrance of a major tributary 

(2) the discharge of an industry or a domestic treatment facility 

(3) an industrial withdrawal 

(4) a dam or other flow-controlling structure. 

The first three categories are necessary in order to reflect changes 
in the quality of water downstream that result from the input of a major 
tributary or an industrial discharge or withdrawal. The fourth category 
is necessary due to the nonavailability of a technique for predicting the 
dissolved oxygen and organic matter below" a dam based on the values of 
these same variables above the dam. A regression method was employed in 
this study using observed values of organic matter and dissolved oxygen 
both above and below the dams; however, the correlation coefficient was 
too small to warrant placing much confidence in the calculated line of 
best fit through the data. As a result, the alternative is that a new 
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reach be created at each dam, with observed values of organic matter- and 
dissolved oxygen input into the model below each dam . 

This inability to simulate the effect of a dam requires that the 
model developed herein re~s tart the simulation process at the ups tream 
end of each impoundment, leaving a gap in the simulated profile at each 
of the dams along the river. 


5.7 Necessary Input for 
Each Type of Reach 

The four criteria outlined in Section 5.6 for designating a new 
reach result in seven different types of reaches that might be encoun- 
tered, with the following distinguishing characteristics: 

(1) SD — ^tream ^Discharging into the river in a polluted condition 

(2) SDH — ^tream ^Discharging into the river and containing _High- 
quality water 

(3) R — Restriction to flow, such as a dam or other flow-control- 
ling structure (this reach is normally given a length of 
0.6 mile) 

(4) N — No input (reach occurs immediately below a dam and is 
designated as a new reach so that BOD and DO values may be 
input at downstream end of type R reach) 

(5) ID — ^j[ndustry discharging into the river ■ (sewage treatment 
facilities also placed in this category) 

(6) IW — Industry Withdrawing water from the river 

(7) IWC — Withdrawing water from the river very Close to 
the discharge point (this reach given a length of zero yet 
still designated as a new reach in order to remain separated 
from the discharge) 

Due to the differing characteristics of the input to each of the 
reaches (withdrawal is considered as a negative input), each type of 
reach will require that different data be input to the model in order 
that the effect of that input might be properly evaluated. Table V 



TABLE V 


INPUT DATA REQUIRED FOR EACH TYl'H OF REACH 
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consists of a summary of the data that would be needed to sufficiently 
characterize the input for each of the seven types of reaches. 

5.8 Division of Warrior 
River into Reaches 

Following the outline set forth in the previous sections for the 
establishment of new reaches, the section of the Warrior River from river 
mile 385.0 to river mile 335.0 was segmented into 25 different reaches 
to be considered in simulating the dissolved oxygen profile of the river. 
Table VI contains a listing of these 25 reaches as well as upstream and 
downstream boundaries and the major input for each reach. 

The actual river mile locations for these reaches were obtained 

from navigation charts of the Warrior River published by the Corps of 
20 

Engineers. It should be pointed out that the numbering convention used 
with the simulation technique employed in this study involves numbering 
reaches at the upstream end of the stretch under consideration and 
proceeding downstream. 

Another convention adopted in this study is, where possible, to 
subdivide the river such that the type R reaches begin 0.5 mile above 
the dam and end 0.1 mile below the dam. The reasoning here is that 
those conditions which would set apart the waters around a dam from 
waters in other sections of the river are: (1) large bottom deposits 

in the forebay, extending perhaps as far as 0.5 mile above the dam, and 
(2) swift tailrace currents, whose velocities might be maintained for 
a distance of 0.1 mile below the dam. 

In some instances, however, the above policy cannot be kept. Such 
is the case around Oliver Lock and Dam, where to extend the reach 0.5 
mile above the dam would include the discharge from the Northport sewage 
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TABLE VI 

REACHES OF THE WARRIOR RIVER USED IN 
SIMULATION OF DISSOLVED OXYGEN 


Upstream Downstream 

Reach R. M. R. M. Input 


lA 

385.00 

382.00 

Ik 

382.00 

367.50 

3A 

367.50 

366.00 

lik 

366.00 

365.40 

5k 

365.40 

361.50 

6k 

361.50 

347.50 

7k 

347.50 

346.90 

8A 

346.90 

346.30 

9A 

346.30 

345.70 

lOA 

345 . 70 

345.70 

lOB 

345.70 

345.20 

IOC 

345.20 

344.40 

llA 

344.40 

344.40 

IIB 

344.40 

344.39 

lie 

344.39 

343.70 

12A 

343.70 

343.50 

13A 

343.50 

343.00 

14A 

343.00 

342.50 

15A 

342.50 

341.60 

16A 

341.60 

338.41 

17A 

338.41 

338.39 

18A 

338.39 

338.10 

19A 

338.10 

337.00 

20A 

337.00 

336.70 

21A 

336.70 

335.00 


Locust &. Mulberry Forks 
Valley Creek 
Yellow Creek 
Bankhead Lock & Dam 
None 

Davis Creek 
Holt Lock & Dam 
None 

Hurricane Creek 

Reichhold Chemical Co. Intake 

Reichhold Discharge 

A. B. C. Discharge 

Intake for Warrior Asphalt Co. 

& Empire Coke Co. 

Warrior Asphalt Discharge 
Empire Coke Co. Discharge 
North River 

Central Foundry Discharge 
Gulf States Intake 
Gulf States Discharge 
Tuscaloosa WTP Discharge 
Northport STP Discharge 
Oliver Lock and Dam 
None 

Hunt Oil Co. Intake 
Hunt Oil Discharge 
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treatment plant, in x^hich case the requirement specifying only one input 
for each reach could not be met. To circumvent this problem, the treat- 
ment plant effluent was included in a reach of length 0.02 mile and the 
reach containing the dam was restructured so as to extend 0.3 mile above 
the dam and 0.1 mile below the dam. This solution provided that each of 
the tX'^o reaches would contain only one input. 

5.9 Additional Equations Used 

In developing a method for simulating the dissolved oxygen profile 
of a stream, the most important equations are those which actually pre- 
dict the concentration of organic matter and the dissolved oxygen 
deficit. HoiN^ever, there are quite a number of other formulae utilized 
in providing input values to the final equations. These supporting 
equations are briefly described in the following sections, 

5.9.1 Velocity Rate Constants 

The three rate constants are comprised of the biodegradation coef- 
ficients kji^ and IC 3 and the reaeration coefficient k 2 • Commonly, the 
value of k-j^, the constant pertaining to deoxygenation of suspended 
organic matter, is obtained as the slope of a semilog plot of organic 
matter remaining versus time. However, in quite a few instances during 
this study, when any BOD measurements at all had been taken, they usually 
consisted only of BOD^ and BOD 2 Q (l-a) . In such cases the value of k^^ 
(base 10) was found by utilizing Equation (10): 

BOD(t) = La (1-10"’^!*') (10) 


which could be modified to the form of Equation (11) . 
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- i log (1- 


BOD(t) s 
La 


( 11 ) 


In those cases where no BOD measurements had been taken by some of the 
industries, the author was forced to rely on some typical values 
given by Eckenfelder^^ for certain types of industrial wastes. 

According to Pyatt^^ the value of venation rate con- 

stant for bottom deposits in reservoirs, can be safely assumed to equal 
0.005. 

Quite a few empirical formulas have been developed for calculating 
k 2 , the reaeration coefficient. Most are of the general form: 


k 


2 


CV^ 

dn 


( 12 ) 


I^Jhere V is the mean velocity, d is the mean depth, and C, m, and n are 
coefficients varying according to the use of each investigator. 

The particular formula providing best results with a certain stream 
varies according to the characteristics peculiar to that stream. Results 
obtained in this study indicated that the equation developed by O’Connor 
and Dobbins^^: 


k2 = 


12.90v°'^ 


d 


1 . 


5 


(13) 


provided the most accurate values for k 2 in the Warrior River. 

5.9.2 Corrections for 
Temperature 

All of the above methods for calculating the velocity rate con- 
stants result in values accurate only at a temperature of 20 ^C and, 
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consequently, these values must be adjusted to the temperature of the 
water: 


ki(T) 


k^(20) 



(14) 


k2(T) = k^cc: 

k3(T) = k3(20) 



(15) 

(16) 


In addition, the value of the ultimate BOD (La) must be corrected 
to the proper temperature : 


La(T) = La(20) 


0.02T + 0.6 


(17) 


In each case. Equations (14) through (17) represent the commonly 
used methods of correcting that particular parameter for temperature. 
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5.9.3 Dissolved Oxygen Saturation 

In order to determine the dissolved oxygen deficit in each reach 
of the river, it was necessary to know the saturation concentration of 
dissolved oxygen for the water in the main stream of the river as well 
as in the inputs of various tributaries and indus trial was tes . 

Quite a few references contain tables showing the saturation values 
of DO at various temperatures; however, in a simulation technique, the 
saturation values are needed in an equation form which is more easily 
adaptable to computer programming. The particular empirical formula 
employed in this investigation for calculating the saturation concen- 
tration of dissolved oxygen in water at a certain temperature was 
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developed by TVA 


23 


and is shown beloxj: 


DOgat = 14.65 - 0.41 T + 0.008 T^, - 0.00008 T^. (18) 


5.9.4 Corrections for 
Changes in Flow 

As each new reach is encountered and the particular input for that 
reach is added to the flow from the reach immediately upstream, it was 
necessary to correct values of dissolved oxygen deficit, and con- 
centration of organic matter for the mixed flow. The dilution method 
is employed in each case: 


^1 Qup ^1 Qin 
" Qup + Qin 


(19) 


D Qup Q Qin 
Qup Qin 


( 20 ) 


La = Qup + .. L a_Q i n ^21) 

Qup Qin 

where Q^p is the flow leaving the reach immediately upstream and 
is the additional flow entering at the head of the particular reach 


under consideration. 



CHAPTER VI 


COLLECTION OF DATA 

6,1 General 

As previously statedj, the primary purpose of this investigation 
was to locate the optimum monitoring sites for DCPs based upon critical 
dissolved oxygen concentration levels. The most promising technique to 
be employed for meeting this objective appeared to be a simulation of 
the dissolved oxygen profile of the river. It became evident that the 
period of the year which would present the lowest 5 and therefore most 
critical^ values for dissolved oxygen would be the summer months of June 
through September, when the volume of flow is lowest and the temperature 
of the water is highest. It also appeared logical that the simulation 
technique and subsequent selection of DCP sites should be based upon the 
critical monthly average values for dissolved oxygen rather than upon 
critical daily values, since the monthly values are less subject to 
variations with respect to location that would normally be encountered 
with daily values. In addition, the data collected was obtained during 
the months of the previous summer, in order to assure that the data 
reflected any recent developments which might cause changes in the water 
quality of the river. 

The nature of a study of this type requires the acquisition of large 
amounts of data for many parameters in order to characterize adequately 
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the quality of the water in a river and in the major inputs. Many 
sources of information were investigated in searching for the data 
needed to develop the model. Information was ob tained . f rom interviews, 
with representatives of the various industries, STORET, the Environmental 

Protection Agency's data acquisition system, and river surveys conducted 

1 L 25 

by Alabama Power Company^^ and Gulf States Paper Corporation. 

It was hoped that STORET, a computer^oriented system devised by EPA 
for ..storage and retrieval of x^rater quality data,^^ would be one of the 
principal sources of information. However, as a result of studying a 
report by Miller and Walters^^ on the use of STORET in water quality 
management and gaining some experience with the mechanics of using the 
system, it was concluded that the data available for the Warrior River 
was not recent, and therefore of little value in providing input to the 
model. 

The remaining sections of this chapter are devoted to a general 
discussion of the sources of information resorted to in the search for 
each type of input data. A detailed listing of the data sources and all 
data gathered for each reach of the river is presented in Appendix II. 

Velocity Rate Constants 

Values of kj^, the overall deoxygenation constant, are necessary for 

each input of industrial x^^aste or discharge from a polluted stream, as 

well as kjL values for the stream flow below each of the dams on the 

river. The k^^^ values immediately below the dams were obtained from BOD^ 

and BOD20 values measured at these locations by Gulf States Paper Corpora- 

25 

tion during “their summer river surveys. In those cases where BOD^ and 
BOD20 values were obtained by some of the industries on the river, such 
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as at Empire Coke Company^ values were obtained through interviews 

with representatives of the particular industry. When such values were 

29 

not known, typical kj values given by Allen and Bodenheimer for paper 
mill waste and by Eckenf elder^^ for other types of waste were assumed to 
apply for the Indus tries in question . 

The reaeration rate constant for each reach was computed by the 

0 ’Connor-Dobbins formula given in Section 5.9.1s while values for the 

13 

other constants — k 3 , r, and s — were all assumed to be 0.005. 

6. 3 Dissolved Oxygen and BOD 
in the River and its 
Tributaries 

The major portion of the data for the parameters was obtained from 

the previously mentioned river quality surveys conducted during the 

2 S 

summer of 1972 by Gulf States Paper Corporation and Alabama Power 

Company. Additional information was obtained from two publications by 

30 31 

the Geological Survey of Alabama pertaining to stream quality. * 

6.4 Industrial Withdrawals 
and Discharges 

The primary source of data on withdrawals from the river for use by 

32 

industries was an interview with Horn followed by a perusal of the dis- 
charge permits on file with the Alabama Water Improvement Commission. 

Data concerning the quantities and characteristics of the wastes 
discharged by the various industries were largely gathered from personal 
contacts with representatives of the industries — Fuller,^^ Davis,^^ 

Woods, Witherspoon^^ and Sandifer^^ — while some additional data was 

obtained from the AWIC files and from historical data tabulated by 
36 


McClure. 
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6.5 Flow Rates 

Principal sources of data regarding average rates of flow in the 

Warrior River were provided by Bowers‘S' and by an impact study covering 

o / 

the effect of the Holt Dam on the water quality of the river. 

Information concerning rates of flow in the tributaries to the 

Warrior River was found in the USGS publication of flow data in Alabama 
o q 

for 1971'^ and in two publications of a similar nature which were 

qi OQ 

released by the Geological Survey of Alabama.'^ ’ Information on dis- 

charge rates from the Tuscaloosa Water Treatment Plant were obtained 
from records kept at the plant. 

^ Water Temperatures in the 

River and its Tributaries , 

Data regarding temperatures of the water in the river were found to 

be plentiful in the previously noted river quality surveys conducted by 

Alabama Power Company and Gulf States Paper Corporation^ during the 

summer of 1972. Similar data for the tributaries v/ere found in Circular 
30 

36, a publication by the Geological Survey of Alabama regarding sur- 
face water quality in Alabama. 

^ ^ Cross-Sectional Areas 

The only recorded measurement of cross-sectional area of the 

qq 

Warrior River was taken above Oliver Lock and Dam by the USGS.'^^ Judging 
from values of x^idth and depth given at this section, it appeared that a 
trapezoidal cross-section could be assumed for the river with the top 
width approximately twice the bottom width at that particular section. 

Since no measurements of cross-sectional area x^ere taken at other 
locations and also as there was a lack of information indicating the 
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general shape of cross-sectional areas, the trapezoidal shape was assumed 
to apply over the entire stretch of the river considered in the study. 
Based on this assumption, the area of each reach was calculated by the 
follox^7ing equations: 

A = d/2(W^ + 1/2 W^) (22) 

which simplifies to : 

A = 3/4 X d X (22a) 

where A is the cross-sectional area, d is the mean depth, and is the 
width of the water surface. 

Values of the width at the water surface at normal pool elevation 
were scaled off navigation charts of the river published by the Corps of 
Engineers . 

Mean depth for each reach was taken as the difference between nor- 
mal pool elevation in the particular reservoir containing that reach and 
the elevation of the river bottom in the reach. Normal pool elevations 
were taken from a pamphlet published by the Corps of Engineers^^ on the 

Black Warrior and Tombigbee rivers. Elevations of the river bottom for 
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each reach were taken from information provided by Bowers. 



CHAPTER VII 


DISCUSSION OF RESULTS 

7 . 1 Comparison of Simulated 
and Observed DO Profiles 

The results of applying the developed simulation technique to the 
Warrior River are shown in Appendix I in a tabular form printed by the 
computer program developed by the author for this model. In order to 
facilitate evaluation of these results, the simulated dissolved oxygen 
values in Tables VII through X of Appendix I were plotted in Figures 5 
through 8 and compared with the profiles based on observed conditions in 
the same stretch of the river. These measured values of dissolved oxy- 
gen were taken from the river surveys conducted by Alabama Power Company 
and Gulf States Paper Corporation.^^ 

Based upon a visual comparison of the simulated and observed pro- 
files, the results, which showed an error of less than 10 percent in 
most places, indicated that the developed model does provide acceptable 
accuracy in the prediction of dissolved oxygen content of the river. 

Some factors contributing to the differences between observed and simu- 
lated profiles were: (1) the highly variable nature of dissolved oxy- 

gen content in a stream, (2) the assumption of complete mixing used by 
Pyatt in developing his dissolved oxygen prediction equations for impound- 
ments, and (3) the. errors inherent in the measured values, which are 
monthly averages of observed values at one point in the cross-section 
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River Miles Above Mobile 

Figure 5. Dissolved Oxygen Profiles of Warrior River (June, 1972) 




River Miles Above Mobile 

Figure 6. Dissolved Oxygen Profiles of Warrior River (July, 1972) 



Dissolved Oxygen at 5-ft. Depth (ppm) 



River Mile Above Mobile 

Figure 7- Dissolved Oxygen Profiles of Warrior River (August, 1972) 
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Dissolved Oxygen at 5~ft. Depth (ppm) 



Figure 8. Dissolved Oxygen Profiles of Warrior River (September, 1972) 
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rather than representations of lateral or vertical profiles within the 
cross-section. 

The lack of observed data at each point was such that the results 
of a rigorous statistical analysis would have been inconclusive. However, 
in view of the fact that observed values varied from the means at each 
point by as much as 10 percent and considering that in most cases the 
simulated values fall x^7ithin plus or minus 10 percent of the observed 
values, the results of the simulation certainly seem to be adequate. 

In those cases where the profiles differed by more than 10 percent, 
the discrepancy can be explained from a physical standpoint. In the 
Bankhead Pool the simulated profiles for each month exhibited the same 
general doxTOx%rard trend (proceeding downstream) shown by the observed 
profiles; hox^ever, the simulated values were somewhat lower than the 
corresponding observed values. This slight discrepancy is considered to 
be due to the complete mixing assumption employed in the DO prediction 
equation by Pyatt, Equation (9), which was applied in the highly-strati- 
fied Bankhead Pool. As a result, the simulated values represent more of 
an average of the higher observed values near the surface and the lower 
observed values at greater depths. 

A sharp difference between simulated and observed profiles is 
noted in the Holt Pool near river mile 362. The observed profiles show 
an almost linear increase in DO between Bankhead Lock and Dam and Holt 
Lock and Dam, x^hile the simulated profiles exhibit a much sharper 
increase in DO immediately below the Bankhead Lock and Dam. This is 
followed by a more moderate increase in DO from river mile 362 to the 
Holt Lock and Dam, The difference between the tx^o profiles is probably 
due to the manner in which the profiles were drawn, that is, by 
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connecting both the measured and calculated DO values on these profiles 
by straight line segments in which the dissolved oxygen values between 
separated points is not elucidated. The higher simulated dissolved 
oxygen value near river mile 362 is x^7hat one would intuitively expect 
as a result of the oxygen-depleted waters released from Bankhead Dam 
undergoing a much more rapid dissolved oxygen recovery than that shown 
in the observed profile due to the large oxygen deficit at this point 
and the corresponding high rate of reaeration. 

In Oliver Pool the two profiles match quite well, with the simu- 
lated values slightly higher in most places. The above situation may 
be attributed to the use of the S treeter-Phelps equation for free-flow- 
ing conditions, Equation (7), which would be expected to yield moder- 
ately higher values of DO in a semi-impounded reservoir such as Oliver 
Pool. 

Excellent results were obtained in that portion of the Warrior Pool 
studied, where the simulated and observed values agreed quite well. 

In general, the model seems to work quite satisfactorily in pre- 
dicting the dissolved oxygen profile of the Warrior River and it is 
anticipated that both this modeling technique and the supplementary com- 
puter program also developed in this study could be successfully applied 
to any river basin. 

7 . 2 Selection of Locations 
Providing Optimum Water 
Quality Data 

Since the primary purpose of the analysis performed in this study 
was to select those points which would provide the optimum water quality 
data with particular regard to the dissolved oxygen content, the simulated 
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and observed DO profiles were utilized in accomplishing this objective. 
The selection of the following points for placement of the anticipated 
five DCPs was considered most advantageous from the standpoint of pro- 
viding useful information: 

(1) R. M. 385.0 (below the confluence of Locust and Mulberry 
Forks) — this point would provide data that could be used to 
evaluate the quality of water entering the section of the 
Warrior River considered in this investigation. Placement 
of a DCP here would provide for collection of data that is 
necessary as the initial input for the system into the 
developed model. 

(2) R. M. 365.0 (one-half mile below Bankhead Dam) — monitoring 
at this location would give an indication of the quality of 
water entering the Holt Pool as well as monitoring the point 
of lov?est dissolved oxygen content in that reservoir. 

(3) R. M. 346.5 (one-half mile below Holt Dam) — Since Oliver Pool 
contains almost all the inputs from industries in the stretch 
of the river studiedj a knowledge of the quality of the water 
entering this area would be of vital importance to effective 
management and control of this highly used water resource. 
Based on this reasoning ^ the placement of a DCP at this par- 
ticular location would permit evaluation of the water quality 
before reaching the highly industrialized sector of the river 
further downstream in the Oliver Pool. 

(4) R. M. 338.2 (forebay of Oliver Dam) — A DCP placed at this 
location V70uld serve a two-fold purpose. First, it would 
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provide for monitoring of the point of lowest DO concentra- 
tion in the Oliver Pool. Second, due to the general downward 
trend of the DO profile from Holt Dam to Oliver Dam, it would 
provide assurance that an improvement in water quality at 
this point would be indicative of improved quality in the 
entire reservoir. 

(5) R. M. 335.0 (approximately three miles below Oliver Dam) — 
Monitoring performed at this location would provide data to 
be used for evaluation of the quality of water leaving the 
investigated area as well as for comparison with data col- 
lected at the entry point to the system for which the model 
was developed. 

A map of the entire length of the river studied showing the selec- 
ted DCP location sites is shown in Figure 9. Detailed maps of the river 
showing the segments used in modeling, the tributaries to the river, 
the locations of all significant waste discharges and the DCP locations 
are presented in Appendix III. 

In addition to the reasons discussed above, these particular loca- 
tions were selected to provide the initial input values of dissolved 
oxygen for each reservoir required in the developed mathematical model. 
Thus, values of dissolved oxygen (and other parameters) obtained from a 
relatively few critical locations may be utilized in the developed model 
to simulate and predict water quality over the entire reach of the river 
selected. As an example, values of dissolved oxygen measured by the 
DCP located just below Holt Lock and Dam (R. M. 346.5) may be used to 
predict DO values in the Oliver Pool as a function of the quantities 
and characteristics of the wastes discharged into this segment of the 



54 


DCP #2 

DCP #3 
DCP #4 


N 


A 



Figure 9. Locations of DCP Monitoring Sites 
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river. Alternately ^ the required DO level at this location may be 
obtained which will produce a desired x^ater quality in the pool under 
varying loading conditions of discharged wastes. In this latter case, 
the information provided x^ould assist in determining the amounts and. 
characteristics of x^7ater to be released through the dam for a given set 
of waste inputs such that the desired x^ater qualic^ x^ill be maintained. 
Additionally 3 the measured DO value and the resulting simulation model 
may be employed to predict the optimum point in time for discharge of 
impounded wastes by the industries in order that the water in the pool 
X‘7ill not be degraded below the desired levels. Similar examples may be 
given for the simulation and management of water conditions in the Holt 
and Bankhead pools based upon the DO measurements at the DCP locations. 
Extension of the modeling technique over the entire reach of the river 
may be performed since x«7ater quality data in an upstream segment serves 
as input conditions to the successive downstream segments. 

Although the model developed in this study V7as based primarily on 
dissolved oxygen, for the reasons previously given, it is in a form 
easily adaptable to simulating concentrations of dissolved constituents 
such as chlorides, sulfates, phosphates, and metals which might be pre- 
sent as a result of drainage from mining operations and runoff from 
agricultural activities. The net effect of many of these materials on 
the pH, conductivity, and temperature of the receiving water will be 
monitored by the DCPs and simulation of these effects may be obtained 
by use of portions of the developed model other than those’ applying to 
the DO prediction equations. For a detailed study of these types of 
constitutents , the DCPs may be moved to appropriate monitoring locations 
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Although not investigated in this study, the use of DCP data in 
conjunction with ERTS imagery may be employed in a basin-wide study 
of the land use-v/ater quality interaction* As indicated previously, 
the DCP x*/ill detect at its location the residual effect of upstream 
surface runoff and drainage from land-use activities. 

In summary, the DCPs are intended tc- provide input data to the model 
developed for the river, enabling the model to be used as a valuable tool 
in formulating and implementing plans for the management of wa.. 3r resources 
and general environmental quality in a river basin. 



CHAPTER VIII 

CONCLUSIONS AND RECOMMENDATIONS 


8.1 Conclusions 

The sites for placement of DCP monitors in the length o^‘ the 
Warrior River from the confluence of the Locust and Mulberry For^vJ to 
downstream of the Oliver Lock and Dam have been selected and proposed 
as desirable locations for the assessment of water quality. Dissolved 
oxygen was selected as the parameter of primary importance and the 
selection of the DCP sites was based on observed and predicted loca- 
tions of critical dissolved oxygen concentrations. 

To implement the selection process 3 a mathematical modeling tech- 
nique and supplementary computer program were developed for simulating 
dissolved oxygen profiles. The modeling technique ^ based on modifica- 
tions of the original Streeter-Phelps Equation by Pyatt, was shown to 
predict values of dissolved oxygen concentration in the Warrior River 
which either agreed quite well with the observed values or could be 
explained on the basis of actual river conditions. Based primarily on 
dissolved oxygen criteria^ the following five locations were recommended 
for placement of the five DCPs designated for use in the Alabama ERTS 
project . 

(1) R. M. 385.0 — below the confluence of Locust and Mulberry Forks 

(2) R. M. 365.0 — one-half mile below Bankhead Lock and Dam 

(3) R. M. 345.5 — one-half mile below Holt Lock and Dam 

(4) R. M. 338.2 — forebay of Oliver Lr-ck and Dam 
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(5) R. M. 335.0 — approximately three miles below Oliver Lock and 

Dam, 

In addition to monitoring dissolved oxygen ^ it is anticipated that 
the DCFs will serve to monitor the other paramieters of pHj, conductivity, 
and temperature to provide information concerning basin-wide land use 
activities. The portability feature of the DCF vzill allow selection of 
other monitoring points for more detailed information concerning specific 
parameters and their relationship to a variety of activities. 

The selection of the specific monitoring sites in this study was 
designed to provide a basin-wide management capability in which the data 
obtained at a relatively few locations could be employed to simulate and 
predict water quality over an extensive river basin area. The selection 
methods proposed in this study were also developed for general applicability 
in the selection of desirable DCF or conventional monitoring sites in any 
water basin. 

8.2. Recommendations 

During the course of this investigation several possibilities for 
further study, primarily concerned with operation of the DCFs after 
installation, have become evident: 

(1) Provide for the collection of supplementary dissolved oxygen 
data at locations where differences were obtained between simulated and 
observed values, permitting modification of the model to obtain better 
agreement with actual observed conditions. 

(2) During the periods betv;een overflights of the ERTS satellite, 
allow the DC? to operate as a conventional automated monitor, transmitting 
data to' a ground receiving station. 
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(3) Include an addition to the DCP apparatus which would allow 
monitoring at different depths of the river, in order to evaluate the 
extent and effects of stratification in the reservoirs. 

(4) Since the DCFs are portable, consider possible relocation of 
one or more of them in order to either expand or reduce the area of 
coverage, such as monitoring the liry of water in a clcv *■: stream 
(Mulberry Fork) in order to obtain base-line aaca on non-polluted con- 
ditions . 

(5) Modify the DCP system for use as a network of conventional 
monitors after the ERTS satellite is removed from orbit. 



APPENDICES 



APPENDIX I 


INFORMATION RELATED TO USE 


OF THE COMPUTER PROGRAM 



Al.l General 


The actual computer programming utilized in the implementation of 

the simulation technique developed during this study was performed in 

/ 0 

the PL/1 programming language. The IBM 360/50 computer installation 
at The University of Alabama was employed in running the program. 

The follox^ing sections contain information directly related to the 
program that would be useful to those interested in the programming 
method and, in addition, examples are given of output for the Warrior 
River during the summer of 1972. 

A1 . 2 Flow Diagram 

Shov;n in Figure 10 is a diagram indicating the general path of logic 
followed by the model in simulating the dissolved oxygen profile in a 
stream. It is intended to show the general nature of the process rather 
than to show the specific calculations made by the program. 

Alo3 Program Output 

The output from the program is printed in a tabular form, with a 
separate table being used for the compilation of data predicted for each 
month. Examples of output containing simulated data for the area of the 
Warrior Basin during the months of June, 1972, through September, 1972, 
are shown in Tables VII through X. 

A1 . 4 Definition of Variables 

Summarized below are definitions for the variables used as input 
data to the program: 
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TABLE VII 


RESULTS OF SIMULATION OF WARRIOR RIVER (JUNE, 1972) 
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TABLE VIII 


RESULTS OF SBIULATION OF WARRIOR RIVER (JULY, 1972) 
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TABLE IX 


RESULTS OF SIMULATION OF WARRIOR RIVER (AUGUST, 1S72) 
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TABLE X 


RESUT.TS OF SIMULA^'LCON OF \JARRI0R RIVER (SEPTEi-iBER, 19 72) 
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lo: 

n 

343o 20 

344o40 

2055o2 


7761 

0*18 

0c.'j3 

0*190 

4*1 

l lA 

IrfC 

3A4o 40 

344c40 

2053. A 






4*0 

l IB 

ID 

344 o 40 

344o 39 

2053.7 


7191 

0*00 

0*201 

0*190 

4*0 

11- 

1 0 

344c39 

343* 70 

2054. a 


719L 

0*15 

0*201 

0cl90 

4*0 

1 2A 

SDH 

343o 70 

343*50 

2 19^. 8 


7462 

0*04 

0*187 

0*200 

4c0 

13A 

ID 

343o 90 

343*00 

2200.2 


7711 

0* 10 

0*186 

0*202 

4*4 

1 AA 

Irf 

343o00 

342* 50 

2172.3 


9072 

0, 13 

0. 186 

0* 192 

4* 5 

15A 

ID 

342o50 

341*60 

2190.8 


7462 

0* 18 

0*187 

0*200 

4*4 

1 6A 

SDH 

34lo60 

338*41 

2 191.^? 


15443 

1.39 

0*186 

0*058 

4*4 

17A 

ID 

338o4l 

33Bo 39 

2192.1 


17589 

OoOl 

0*186 

0*030 

3*2 

1 34 

R 

338*39 

338*10 

2192. 1 

2320c0 

289L2 

0*25 

0*128 

O 0 O 22 

3*2 

194 

N 

338o 10 

337*00 

2320.0 


9562 

0*28 

0*065 

0*130 

6*7 

204 

U 

337oOO 

336*70 

2319.3 


7779 

0*05 

0*065 

Oc 140 

606 

?1A 

ID 

336o70 

335*00 

2320.0 


6693 

0o30 

0o065 

0*150 

606 


00 

cr:D 


5.;? 

5o0 

5o7 

3o0 

5o0 

^o 2 
<^c2 
^.2 
^^o 1 

4c0 
4c 0 
4«0 
4o0 
4o4 
4« 5 
4c 4 
4o 4 
3o2 
3o2 
6o 7 

6o6 
6o 6 

6c5 
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TOT_REACHES 

total number of reaches into which the river is 
divided » 

NOR 

num.ber of reservoirs on the stretch of the river 
being modeled. 

NOM 

number of months for which simulation is to be 
performed. 

K320_E 

deoxygenation constant due to bottom deposits at 
20 . (base e ) 3 days”^. 

X 

deoxygenation error term, days""^ • 

Y 

reaeration error term, days~^ . 

LAI(S) 

concentration of organic matter in reservoir S 
before waste inputs from upstream are added, mg /1 

DEFI (S) 

concentration of dissolved oxygen in reservoir S 
before waste inputs from upstream are added, mg /1 

FORMCS) 

type of equations used for predicting DO in reser- 
voir S, Streeter-Phelps or modification by Pyatt 
for reservoirs. 

TYPE 

category of input for the reach in question — SD, 
SDH, R, N, ID, IW, IWC. 

REACH 

number of the reach. 

MONTH 

month for which the data was taken. 

UP STM 

upstream boundary of the reach, river miles. 

DrasTM 

downstream boundary of the reach, river miles. 

DISCH 

discharge of a tributary, cfs. 

DIS_I 

discharge of an industry, mgd. 

WTHDR 

amount of withdrawal by an industry, mgd. 

REL 

amount of flow released from a dam, cfs. 

RIVTE14P 

temperature of the river water, ^C. 

INTEMP 

temperature of water in the input to a reach, , 

WIDTH 

surface width of the river, ft. 

DEPTH 

mean depth of the river in a reach, ft. 
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30D5 

BOD 3 of the input waste at 20 mg/1. 

B0D_pV7N 

BOD 5 of the river water (20 ^C.) immediately 
below a dam, mg /I. 

DO 

dissolved oxygen concentration in the input to 
a reach, mg/ 1 . 

DO_DI'JN 

dissolved oxygen concentration of the river water 
immediately below a dam, mg/ 1 . 

K120JL0 

deoxygenation constant at 20 . for the river 

water or input to a reach (base 10) , days 


A1.5 Listing of the Program 
Statements 

On the following pages the actual PL/1 statements of the program are 
shown. It is hoped that any potential user employing the information 
provided in the preceding sections of Appendix Ij in addition to these 
program statements, could either utilize the program intact or easily 
adapt the programming to fit his particular needs. 



PROC OP T I Ofv! S < M A i N ) ? 

/ V« <s ☆ »it >:« **^ * *4; :;{ V V si s*,e st: C; vJ O ^ ^ vt O V ❖ 

/ *:i: >;< I,’: V i\z >|c 4 >»^ >:«>;: i|4 i|r k;? :ts 4- * 4.' sU 4< * i/: 5^: 4: 4;: >>: j> i^: s;; <:<: <t sjr^j jjc jji C: <£ ^tr^O^nCrOCrs^ :^:4:4i 4 j;;j 0 <* v >^iCc J> / 

THIS PROGRAN*: S;'':ULAYES DISSOLVED OXYGEN IH A RIVER ENVIRCNr.ENT 

/ *>:c -4: s;; >;i <: ijc 4: << j;c *4's;; -4 *4: 4s ‘4 4 : iC; 4: <e vijf 4: * 3i< * * 4: ❖ 4; <; 4i 4: ❖ 4= v£ ❖ 4 : <2 * 4: 4:4* ❖ 42 4*4:4*4* / 

/>js<tx;<i^^»t^&<si5s)(4i:sj‘^>:s4:4:*C24:^4<4s4:4:4<4f4:4:4:4:i>:5jE4«5y*4*5i*4<4:4‘5Ct4:4f<:’v:4:^<*4*4s4*4c>;s4:4:4:v:C*4*4:*<* 4** 4:4* C: 4* 4* ☆ <:/ 
nCL ZILCH CHARJ L3 S 

OCL ^EXPI,EXP2) FIXED DSC I ^UL ^ 5 9 3 ) ; 

DCL FORS{^^; CKARIG: VARYING; 

OCL CMT FIXED DEC I ?->AL ( 3 5' ; 

DCL CAPA:iTY(i65 FIXED DEC I KAL I 5 9 0 ) ; 

OCL ( 16) yO^OUTI 16) ) FIXED DEC 1 MAL { 6 9 I ) ; 

OCL {X^NSWCHiO FIXED 0ECIMALC 5 9 3 ) ; 

DCL (TYPEj^FACH) CKARI3) VARYING; 

DCL ■: MOJ 9 MONTH) CHARI 9: VARYING; 

DCl DATE(^^: CKAR( 15; VARYING; 

DCL (rOT_.REACKES 9 NOM 9 NCR 9 HvN 9 S 9 T) FIXED DEC I MAL ([ 2 I> ; 

DCL ( SLOPE 9 intercept; FIXED DEC I HAL il 6 3 ; 

DCL {CVD 9 E 9 F 9 G 9 H 9 I 9 J 9 X 9 L 9 O 9 P 9 Q) FIXED OECIMAUTp^jD ; 

DCl ( 500_UP 9BOD5 9Lk_I N9LAyL_DUN9LA_IN:T9 B005_I NIT sBOD^Ok-^NvLAZOo 
LAH<^>5 FIXED DEC I MAL J 6 9 1 5 ; 

DC. { R IVTEMP 9 I NTESP9D09DO_UP9 0G_DIJN9 0GSATR9DOSAT'J9DEF_UP9 0£F_DWN9 
OFF^INjDEFv DEF^INiTvOCj^INI TffDEFIi^,^ n FIXED DEC I LUL I ^ 9 i D ; 

DCL ( KI20_UP9iCl20_10 9{<120_E9K220_lO 9K220_E 9K120_HIX9Kl9l<29{<39 
K320^E9W?Z9A9S; fixed DEC! MALC 5,3) ; 

DCL ( Q^UP 9DISCK9l*IDTH9DEPTH9FL0UvREU FIXED DEC IMALfi 69 1 3 5 
DCL { UPSTM 9 DWNSTMyOI S_l 9 WTKDR ) FIXED DEC I MALI 6 92 5 ; 

OCL (LENGTH 9 CSA: FIXED DECIMALI 6 ); 

DCL VEL FIXED OE C I M AL ( 9 2 5 ; 

OCL TIME FIXED DECIHAU 592); 

DC) VOLUME FIXED DECIMAL! 105; 

M0{ 1 ) = s JUNE = ; 

MO( 2) =® JULY^ ; 

M 0 ( 3) = '^ AUGUST® ; 

M0( a; =’ SEPTEMBER® ; 

GET lDIT! T0T_RE ACHES 9 NOR 9 NOM9!<320_E9X9Y9Z ILCH5 (C0L( I 5 9 F I 2 I 9 COL <l 2 U 9 
F! 2 ) vCDLiOL; 9 F i 2 ) 9 COLi^;^i 5 9 FI 593 S 9 COLI 5 U 9 F {[5935 9 COL 2 6U 9 Fd 5 9 3 5 9 

X ( 1 ) 9 A 2 L ; 5 ; 

00 S = 1 TO NOR; 

GFT FDmLAUSUDEFHS; 9 FORM! S) 9 ZILCH) ( COL I U 9 F 2 5 9 1 ) 9 COL 2 10 ) 9 F 2 ^9 lU 
C0LC20) 9 A( 3 > 9 X 2575 9A( n ) ; 

FND; 

M=I ; 

^FSET; Q_UP = 0 . 0 ; 

X120_UP=0o0; 

0£F_UP=0o0; 

ROO_up = o oO; 

PUT PAGE; 

P'JT SKIP! 13) EDIT :« SUMMARY FOR VJARRIOR RIVER®) 2COL2^8)9A5; 

PUT SKIP22) EOITC®Ro Mo 335o0 TO Ro 335«0®5 { COL 2 <^7 ) 9 A ) ; 

DATEr;: =MC( M) j; 3 9 1972 ®; 

PUT S<IP22) EDITCDATEIM) ; 2 COL 2 54^ 5 9 A ) ; 

PUT S'<IP(3) E0IT2®REACH®9 ®TYPE ® 9 ® UPSTREAM ° 9 ®OQUNSTREAM® 9 ® FLOW® 9 

®RELEASE®9®CSA°9 ® travel ® 9 ®K1 ® 9 ®K2 ® 9 ®00® 9 ® on® ) 2 COL 2D9A9COL2859A9 
CDL2 155,A9C0L{25) 9A9C0L (38) 9A9COL2^6 ) 9AcCOL2 59) 9 A 9 COL 2 6 8 2 9 A 9 
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cni ( 80 : 5 COL ; 91 ) 9 A,CO;. { 102) V A vCOL(114 5 9 A5 ; 

^uT SKIP EDI n ' NG« « V ® Ro n;<.%’Ro s m CPS ) « v ' { cfs ) % m sQo fTo $0 ti:eie% 

MGASE 6 )S’iBASE E 5 % 'UPPER E NO ' 9 a OirJ E R ENO^ ) ( CGL i 2 ) 9 A ^COL i 1 7 ) A 9 

CDL{2B)9Av:nLC3839AvC0LK75 9 A 9 COU 56 5 9 A ,001(69) ,A,COL ( 77) , A, 

CGU 38 ) 7 A, :0L( 99) V AvCOL I 1 1 , A) ; 

PUT SKIP E!)ITM { DAYS) « 9 M PP^O S M PPW) M ( COL ( 68 5 , A , COL { 10 1 ) , A , 

CDL( I 13 ) 9 A ) ; 

P‘JT SKIP! 2) EOtTC «) (C0L(2),AJ; 

NJ = 0 ; 

T = 0; 

AOO: N = NJ-}'I ; 

if N = i THEM DC; 

T = T V 1 ; 

LA^ INI 7 = L AM T5 ^ 

OEf^INIT = OEPi:m 
END ; 

ELSE on TO check; 

CriECC IP N > TGT^REACKES THEN DO; 
if S=NDS THEN GO TD LAST; 

ELSE GO TO NEXT; 

NEXT: N-Mvl; 

GO TO RESET; END; 

ELSE 00; 

READ: GET E 0 1 T ( T YP E , R E ACK , HONTH ) {COL( U , A( 3) , COL ( 6 K A ( 3 ) , COL U L P 
AI9) ) ; END; 

/ V VC jjj sj: <: >;c * *4; ^ ^ :1: i;« ☆ V *4^ ❖ V 3\i »/: JC: 5*£ Jji >> ❖ 3^: ❖ * 3jc Ci jj: :CJ * ❖ 

/ >;£ Tjc 4c 4c 4t si: i;c 4 : <c VC v 4c V 4c 4: C: 4-’ 4c sCc s:c 4: 4< 4 : 4: 4c s); i): « <c 41 ^c <: 4 ; jJ: j>: s;: i>: ^ 4c 4: 4; 4< 4: 4: 4; V- 4= fr O <c 4 c / 

/4c>;c4:4:4c*;i<c4c4c4^>fr4c4c* POLLUTED STREAK DISCHARGING INTO RIVER Cz4c4c4:4:4<4!<;i:csCcCc4:04:/ 

/ 4c 4s 4 : :^c 4c 4c 4: 4c lit 4c * 4c 4 : 4 ; 4« V V '4c ❖ 4c 4: 4 : sC: >;c 4« 4i ij: <c 4c 4c j(i sis 4e sfc 4s sj: <: 4 : 4 : 4 : 4s *5: 4: 4 c 4: 4c 4c sCc =Jc 4c ❖ 4c 4c 4c 4n;c Cc O c^c 4: 4c / 

/ V 4^' 4c vC 4c 4c 4* sjc V *4c 4c 4* 4c 4t 4 * V*' V '<c \: 4 c •,< s^ V* ^c 

IF VVPE=»Sn« THEN 00; 

ST^ EAS: GFT E0 1T( JP STM , OHNSTH , D I SCH , RI VTENP » H I OTH , DEPTH 3 tC0Ll25)sF26s2Jt 
C3L( 3S; , r( 6,2) vCOL H-55 vEI6t 1) sCOLt 57) , F ( 4 , I » , COL ( 65 ) »F(6» 13 , 

COL( 75:,F(6,i3 ); 

GET EOIT( 51D5,00,!C120_IO,INTEHP,Z:lCH) IC0U25J , F ( 6 , 1 ) , C OL ( 37 5 o F J <>» II , 

C0L(66i ,F( 5,33,CnLii57i ,F{69l) ,X( 19),Ami ; 

IF N=1 THEN 00_UP=00; 

ELSE D3_U“=0nSATR - DEF_UP; 

OOSATW = 16.65 - C <> 6 1 r NT E HP v 0 « OOB* I NTE MP>i*2 - 0. 00008* I NTEKP<'P3 5 
_ENGTK = (JPSTH - OUNSTHS * 5280.0; 

CSA = 0.75 <= v;:dt;-; * depth; 

VOLUME = ;Sfl V length; 

FLOv) = 0_UP <- DISCH; 

VEL = flow/ssa; 

TINE = a.ENGTH/VELi / 36600.0; 

XIZO.E = K120_10 * 2o303; 

KI20_SIX = ( ( < 120_U?*g_UP ) ^ { K120_E*DISCH) ) / FLOW; 

K220_10 = (12.90 * VEL**0.5I / ( OE PTH** 1 . 5 ) ; 

K220_E = K220_10 * 2.303; 

<1 = K12 0_HIX * 1.067**(RI VTEHP-20.0) ; 

<2 = K220_F * 1.0159**( RI VTEHP-20.0) ; 

<3 = K320_F * l.067**(R! VTEHP-20.0I ; 

LA_IN = B0D5 / C 1.0-1 .0/; 2.71**! 5*IU20_E ))) ; 

La20 = ((BOD_UP * Q_UP) •> !LA_IN * 0 1 SCH ) ) /FLOW ; 

LA = LA20 * (0.02 * RIVTEHP <■ 0.60J; 
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DEF.-rv) = onsATU - do; 

03SATR = - 0.<ri*RI VV£'-‘? -o- 0 . 008*R I VTEMP‘>’>2 - 0.00008*81 VTE^PC“>3; 

C£F = U ( nOSATA-DO_U? ; * 0_UPS ■«■ £OEF_IN * OISCHS) / FLOW; 

:F FjR'‘-^(T* = "RES' THEN 00; 

= FLOW « (86^>oo.o/vqlu;-.e: ; 

IF N= . THEN Z = OlSCrl * i «,6<-00. 0 /VOLUME ) ; 

FLSE Z = Q_UP * { 86<>00.0/V0LUnE ) ; 

A - K 1 V i< 3 W ; 

B = R 2 ■«- w; 

B30_S0: L_OWN = L A_I NI T* II . 0/ { 2 . 71 *« { A*T I ME H J J <>■ {IZ*LAI/A!I * 

C loO-l oO/( 2.71**{ A*Ti:‘.£ 5 ) 5 ; 

C = D6F_INrT - IZ*r-cF5/B: 

0 = / ( A-B5 : 

6 = LA_iNIT - 'Z*LA;/8; 

F = 1 oO/{ 2.71**C B*TIMc 5 } ; 

G = KL/{ B-A ) ; 

K = LA.INIT - {Z*LAJ/A; 

1 = l,.0/{ 2.71**£ A*TIMEn : 

J = 2/5; 

< = OEF V (Kl*LA)/A; 

OEF_SO: 0FP_3V.'N = £C-e-£DvEn*F (J*K5; 

end; 

ELSE 00; 

L_5i.‘N = (LA - ;:/:U) * c 1 oO/( 2.71**(K1*TIME) >) <• X/Kl; 

CXPL = / < 2. ?!**{ K1*TI me; S ; 

EXP2 = lo'!' / ( 2.7i'-;'*( K2*TIME} 5 ; 

IF R2 = Ki THFN 00; 

L = Kl' « TIME ; 

1 = (LA-X)/Kl; 
o = (XvY)/Kl; 

OEr_Dh'N = (L^n -c- OEF -fr P 5 * EXPl - P; END; 

ELSE 00; 

0£F_0V1N = C(K1*LA - X5 / IK2-i<UJ * (EXPl - EXP2K > UX«-V)/t(2) * 

£ I .0-1 oO/i 2. K2*TIKE ) > J -c- OEF* ( I .0/ ( 2 . 7 1 ** ( K2*T I HE } D J ; END; 

END ; 

IF LENGTH < ( 0 . I * 5? 80 . 0 I THEN 00_0WN = DO_UP; 

ELSE 03_DWN = OOSATR - 0EF_01.'N; 

PRINT SO: P'JT SKI? EO I T ( REACH , TYPE ,UPSTMv OUNSTM , FLOUo CSA» T I ME ffiCl<,l<2 , 
00luP.70_0WNI JCOLJ 2 S V Ai35 ,C0L(9J , A ( 3 5 , COL (! 1 7 } , F t 6 , 2 J , COL { 28 5 o 
F{6,23«COL(37},FJ6vlj ,COL{ 583 ,F(6) ,COL(68 5 vF(5,2l » COL ( 79 ) , F ( 5 , 3 J , 

COL (90! ,FI 5,3),C0L( 101> , F £ v 1 1 ,COL ( 1 1 3 ) <, F , 1 ) J ; 

3_JP = flow; 

BOO.JP = L_0WN / (0.02 « RIVTEMP 0-6); 

OEF_UP = OnSATR - 00_0UN; 

K120_UP = ki20_m:x; 

GO TO AOO; 
end; 

y 5^ Gjf VI 5jj 3^ ^ V' V ^ V C* W V *.1 V iji ;}c VI Vr 41 J}< ^C* V ^ ^ ^ V v* lil v Cl Cl C* C* V 7^1 Cl C* Ci **i 

/ VC Cl C: l: i;: C: C« C- Ct ❖ 1 : 1 * C*C: ❖ <c C^ 0 * C? C* C:* « Ct dCJCrOCi 00^0 *0 * C: C: C: Cr <: <: / 

/,*,ci:x3;c;;:.;:V-C:C:<ii;iC:CrCi'V GLEAM STREAM DISCHARGING INTO RIVER 

/ aj: V CiCi Ci v' *C: C: *>%' d i;; V C« V V C^Ci Ci^C? CJ VOd ❖ C: ❖ CcC** Cr « * ^CjC: ☆ C: C^ CtC« C^Ci C* C:* C: * ❖ Jj: * CsC: C«/ 

/ :;;v C<*C:CiCiC:A>!c sic aUCcj^cv* ^cCtCtCiCiOCtCrOCiCtCiCtCiCtCtCiCcCs Cl/ 

ELSF IF TYPE = =SO:-i« THEN DC; 

STREAM CLEAN: GET EO I T{ UPSTH, OWNSTH, 01 SCH , R I VTEHP , WI 0TH,O£PTH 5 
(COL( 25) ,F( 6» 2 5,C0L( 35 ) » F £ 6 , 2) , COL ( 45 ) ,F ( 6 , 1 K COL ( 57J ,FI4,IK 
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CC*_ vo^rf.'*\ 05 ±»oCO»_'i 5 

GET- EOKiDr). :.N!TEn;? 5Z”LCKJ {001(37} ,F(<>,1} ,C0L{57} »F(l<>»llil jKl 19} jAtl 3 5 5 
IF \'=: THEN On_UP = DO; 

ELSF 00_OP = 0GS4TR - DE;-_UP; 

D3S£.T'.; = - Oo-fjl^INVEHP «■ 0, 008* I NTEKP«==?2 - 0»00008*R1VTEMP**3; 

LENGTH = {UPSTG - Dl.'MSTr-",} * 5280o0; 

CSA = 0o75 * '..'iDv;-; * depth; 

VOlU' iE — C SA * LEi'JGTK; 

FL?)'..' = Q_GP <- OISCK; 

\/EL = FLO'v.' / CG;,; 

TINE = (LENGTK/VE-3 / G6<-00.0; 

K120_SIX = :XI20_UP Q_UP i / FLOVJI 
K220_10 = 112.90 * VE1**0.55 / ( DEPTH**1 ,5 } ; 

!<220_E = K220_10 * 2o305; 

Kl = l<120_!-'.!X V 1 oO‘'> 7**:R! VTEr-'<P-20.0) ; 

;<2 = ;(.220_E * loC;.-J9**;RI VTE.NP-20.0) ; 

KZ = K320_E * loC47**IRIVTEr-:P-20<,03 ; 

LA20 = iBCO_GP * Q_U?S / FLOW; 

LA = ^A20 * v0,02 * RIVTEHP o 0.605; 

GEF_I.y = GOSATl.’ - DO; 

DOSaTR = K-.65 - 0.6i*RI VTE;;P 0.008*RI VTEEiP**2 - O.OOOOS*RI VTEKP**3 
OFF = {{{OOSATR - D0_'JP3 Q_UP5 <- (OEF_IN * 0ISCH3J / FLOW; 

IF f3rn;t:' = «res^ then go; 

i; = FLOW * : GoNOOoO/voluhe; ; 

IF N=1 THEN I = DISCH * J 86400. 0 /VOL UHE 5 ; 

ELSE Z = 0_UP * (G6400o0/V0LUME3 ; 

A = <I V K3 w; 

8 = K2 W; 

800_SOH: L_DWN = :.A_I HI T*{ 1.0/J2. 71**« A*TIKE5 } } «• {«Z*LA)/A} * 

( l.0-L.0/{ 2o71**: A^TIHE ) ) ) ; 

C = DEF_INIT - {Z*DE.-5/B; 

0 — Ki / ^A S/o 

E = LA_INIT - IZvLAJ/G; 

F = loO/{ 2.7i**(B*TIKE3 5 ; 

G = ICl/IB-AI; 

K = LA_INIT - M*LA5/A; 

1 = i o 0/ { 2 . 7 1 { A* T I HE 3 3 ; 

J = Z/B; 

IC = OFF (KI*L*^5/A; 

OEF_SOH: OEF_CV.'N - { Co-C Ov£ S }*F -a- IG*M3«I «• {J*X3; 

END ; 

ELSE 00; 

L_Ol.'N = ,'LA - X/K15 * I 1 . 0/ I 2 . 71 ** ( !<1 *7 IME 3 } } K/Kl; 

e;;p: = i.o / i ; ; 

EXP2 = 1.0 / I 2.71**I!<2*TrM£} ) ; 

IF :<2 = Kl THEN GO; 

L = XI * TINE; 

0 = ILA-X)/X1; 
p = Ca-£-y)/ki; 

0£F_3WN = (L*0 V OEF i- P3 * EXPl - P; END; 

ELSE DO; 

0EF_0'.-.'N = (CK1*LA - X5 / i:C2-Kl!) * (EXPl - EXP25 v ((X<-V3/K2} * 

(1.0-1. 0/(2. 71**(X2*T:;-‘,En} <- 0EF*(l.0/(2,71**(K2*TIf;E3 ) SS END; 
END ; 

IF LENGTH < (0.1*5280.05 THEN 00_0WN = 00_UP; 
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ELSS DO_DV.’M = Or^SATR - DH?_D'.:W; 

PRINjr_SDK: PUT SK'iP S 0 1 T ; RE ACK , TV PE , UPS TM 5 D'.JNSTM 9 FLOli o CS A 9 T I ME «. K 1 5 K 28 
DO_JP, Dn_DiC!\:; (COL^E:,;. 'SisCOLCg! vACSD sCOLHITII 9F(I692!)»C0L£281)9 
F( 6 , 2 ) fCnU 3T; ,F( 69 IS 9 COL: 585 9F(65 ,C0Li68S ofi5,2Si ,COL J 79 ) , F4 S, 3 J 9 
C3L( 98 5 ,FJ 5,3 } 9 COLI: lOi 5 , F 4 4 , 1 ) , COL U 13 } , F K , U J 5 
Q_U? = FLOU; 

p.3D_UP = L_D',;:; / ;o.,02 RiyTE.'-5P v o<,6t; 

OEF_JP = DOS.'.TR - DO_D'.,'W; 

K120_UP - K120_M:X; 

GO TO AGO; 
c\'D ; 

/ V C >:•• V V <: ^ !> O * O V V O ■.-• ••:; O-U * ■'/ O * O * C= * * V C; C-- * * * C: s: 15 j:; * O * * O O * / 

V ^ V ‘fi ^ V* C' V ^ ^ ^ O 't* ^ ^ V* 4 ^ V '<• O V V V ^ N* V C* ^ ^ ^ V O / 

/**-<* reach C 0 ;'JTAIHJIVG a dam or other flow controlling structure 

/ \' VT V ly* ^ i'll i,i i,J 1,1 i.i i,i i,i i|i i»i i|i V '.t ill i|j v* V- *r* V V ^ ^ v- ^ V V V V ^ C* 4* V«- C» ^ 

/ V j;: V V ilJ V :|J v ijr C: ij: il: sC; ;;. sji l;i C: jCt V* 5:^ V * * ^ * s;? * * <•" C* / 

ELSE If TVPE= ^R- VXETC DOv 

RESERVE I R: GET E 0 I 7 CUPS VM v GL'N STf-^ 9 REL RI VTEMP 9 W : DTKo OE ?TK i iCOL (125 5 9 

F{ o, 2; , COLi 355, F « 6, 2 5, COL I 45 5 , F {6,1), COL 5 57) ,F44, O ,COL 165) 9FC69 U 9 
COL( 75; ,F: 69 15 ) ; 

GFT E D i T i B no_Dl,'N , OO.GUN , !< 1 20_1 0 , Z I LCK 5 { C OL I 25 ) 9 F I 6 9 U 9 COL I 37 ) 9 F (! 4 9 1 ) 9 
COL{ 56) ,Fv 5, j / , 55 1 19) vA( 1) ) 5 
IF N=1 THEN 00_UP = 00; 

ELSE OQ_UP = D0SA7R - DEF_UP| 

LENGTH = (UPST.4 - DUNSTI-1 5 -7 5280,0; 

CSA = 0,7 5 ••■' WIDTH v DEPTH; 

VOL USE = CSA V LENGTH; 

FLriW = 0_UP; 

VEL = FLO'.-i /CSA; 

TISE = (LENGTH /VEL) / 86400,0; 

<i20_E = K120_10 V 2o303; 

!<;20_Kix = ;<i20_ir; 

K220Z10 = (12,90 p VELPvOoG) / «D£PTH4*l<.5) ; 

K220_E = K220_i0 v- 2,303; 

<i = !<:20_,-;:x 4 i„047pp;r; vtemp-20„05 ; 

K2 = K220_E l,o:59P>^;r<IVTEHP-20oO) ; 

K3 = .<320_F P 1«047*--:-Cr;vV£NP-20,0) ; 

LAZO = 800_UP; 

LA = LA?0 4 (0,02 4 RiVTEr-5P v 0,60); 

DOSATR = 14,65 - 0, 414R I VTEF5P <>• 0 , 008 4R ! VTEHP4*2 - 0, 000084RI VTEMP443 ; 
OFF = DOSATR - 30_UP; 

DEF_Ol.'N = DOSATR - DQ_D'.:N; 

PRINT_R:PUT SNIP FDIT :rEACH 9 TYPE, UPSTH, DWNSTMs FLQH 9 REL 9 CSA 9 TINE 
, i<l9 R2, 00_UP, 00_0UN) uC0L{2),A{3), CQUO), AC3), C0LJ17J9 FJ692) 
9 COLCZSS, F( 69 25 v C0L(37)9 FC6,1)9 C0LI475 9 F«69l)9 C0H58)9 F{6)9 
C0LI68), F( 5,2 5 9 C0L(79)9 F(5,3)9 COLiOOjj F { 5 9 3 ) 9 COL { 101 ) 9 F(4ol»9 

CtJLI ,13)9F(49 ji) ) ; 

0_Uo=REL; 

DEF_UP=OEF_OHN; 

Ki20_UP=K120_SIX; 

GO t5 ado; 

END ; 

^ sjj 3|r ij* iC^ V.I j|i S|j s/ V V V ',* V 4* V 4* ^ i{j l|l s|J i|i Sy! iyi !*, iji 

f *i;s ;:;4: V Ct i}( V v v|: s|{ ili ',11: i/i ^ O i: j,;;:; * * s|: * :C! jC; * :<: * j;: Cj jI: j;: <£ :Jc <X jC: << 

REACI-; CONTAI i\JI NG NO UNPUT TO RIVER 
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-4, qccuRS EELOW REACH CCNTAiWIl\:G A 0AM 

/ 4c »;•>;: 4: ^ i|f 4e * >!t * 4' V 4; v'4c 4= O 4^4; 4c 4-' 4: * 4c 4i * 4: * 4; :>4: O si 4: 4* v 4c ❖ <!s;i v 4t ☆ 4s 4t<i 4* 4s <i 4i C: O <5 0 / 

/vOOO^Vsc0 4s'.,c4ci|c4i4cc<c4c4c4cc,ccyCi,cc,c4c‘.^;{ii^\;4zc^i^jjc4i<i4c-4c4sc4i4ccicc^4c4s4scCcvcJcsiccCsslccIc4i4s4ssCi*»sOO'vccii<sc»sOcSc4c4s4s4i4c4i/ 

else IE TYPE = ‘i\!= TKEiS DO; 

T = V -;• i; 
l.A_INn' = LAI( t; ; 
nEF_!M r: = oee i ; t; ; 

M0_c'>j7: GET EDIT <GP5TM„ OL'NSTMo DISCHs, RIVIEMP, WIOTK, DEPTH! JCOLt 
251, F(652>9 ZCIIZ-jU fib,!), C0LU5D, F565lSo C0ILt57}» F(<>,11)(, COL { 

65/9 Fi09L/9 ^GL /p)/ t"v69i.))9 

GrT F0tT;K120_:0,Z-\.C:-:;CCCLJ26>, F(5,3I 9XU9} pAUS 5 5 
cF Mssl THEN Dn_UP = GG°o 
ELSE On_UP=sDOSATR-DE.-_i;P; 

LFN!GT;-I=C UPSTK-Dl.’MSTM', 05280.0; 

CSA = n„75 o '.-hot;-; >:s depth; 

VOL USE = CSA * LFiVGTH; 

FLnK'=sO_GP-:-D:SCK; 

VFL=?Ln'.!/CSA; 

TI YF= ; LEtvlGTK/VEL 3 /CbbCO^Q', 

<120_E = :LI2C_1002. 305; 

K 1 2 0 _H I " = K _ 20 _E ; 

<220_10 = ': 12.90vVEL0c:c0o5S /CDEPTHOoi ,5} ; 

'<220_E = K2 2 0_10O2., 305; 

K 1 = K ' 2 0 _ S I * 1 o 0 <9 7 * * ; R i vy Y E r--, p - 2 0 o 0 5 ; 

<2 = ;C2?0_EO 1.0150>;sc-(;r j v/Yt:',P-20.G5 ; 

<3 cr K320_E o 1.0<j7o*:c;ri vTEHP-20.0) ; 

LAZO = B0G_D1JM ! { 1 . 0- I . 0/ I 2 . 7 1 * o ( 50|< 120_E i J) ; 

LA = LAZO 0 (0.02 0 revte;-'.? ■> 0.60); 

OOSATR = L<9.65 - 0 « <;■ I OR I L'TEHP + 0 , 0080R I VTE MP<s*2 - 0. OOOOBORI VTEHP-^03 ; 
DeF=OOSA7R-DO_UP; 

IF pCR*-',(T; = »RFS“ THEM DO; 

W = FLOi.' o C S6 c!/00o0/V0LL;KEJ ; 

Z = 0_|JP * ( S6400.0 /VOLUME ; ; 

A = K 1 ■;■ i< 3 -s- w ; 

B = i<2 + w; 

BOO_N; L_0'.;N = L.c_;HITO{ :.0/(2.71*0{AOTIME1) ) J UZ*!LA)/A) «= 

i 1 oO-i .0/1 z.tlooImOtihe} 5 J ; 

C = DEF_Im;t - :Z0DEFI/3; 

D = <1 / ( a-b; ; 

F = LA_IN!!T - I 2'OLA I /D ; 

F = i.0/( 2.71 *o’£*VIHE’; 'J ; 

G = R1/(B-a;; 

H = LA_IMTT — CZs.cLAI/A; 

I c= : „0/( 2.71YOJ ./VT'iHE J ) ; 

j = z/b; 

K = OFF 4- {K10LA5/AS 

OEF_N: DFF_DV.'NO =c= I C-:- 1 D0£ I I *F ■>■ IGOHJYI 4- JJOKJ; 

FMD ; 

cL S tr DO ; 

L_D'..’M = ILA - R/ICIS O Ci^O/IZ.yiooiKlOTJME} )! 4!- X/Kl; 

ERPl = 1.0 / ( 2.7 :o--;si;;li<cv:k£) 5 ; 

FXP2 = 1.0 / ;2,71ooi;L2<cT:i-'.E) J ; 

IF K2 = Id THEM CO; 

L = Id p time; 

0 = ILA-K)/Kl; 



t’ = (x<-yj/!<l; 

0f<^_0..'N = (L*Q V DSF ■> P} =•■= EXPl - P% END; 

EL'^E ot; 

n5F_0',;M = - X) / ::<2-Ki;f * JEXPI - EXP2) <• ( lK-:-Yl)/!<2) * 

( ioO- UO/( 2«71 -x*{k 2--^T; '--.Ei ) 3 <- DEF=> J 1 »0/ H 2.71<“X J K2«‘T HME J J J 5 ENDS 
END : 

!F LEN"jTS < (0.1<^5280<.0; THEN 00_DWN = DO_UP; 

El. SF 0.1_OV-3.N = DnS/-.VR - DE.-_DWN; 

PRI NT_N:PUT SiCl'P EDI V : R EACH » TY PE , UPSTt-'.j OV3NS VMc- FLOW , CSAs Til KE sKl s X2 9 DQ_i;? 
,o5_7)CN } CCOL I 2) vAIS) ,C0L? <?’/ 0 ACaUCOL ( 17) ,F(6j2nC0Lt 28)oF!6,2) 9 
COL ( ^7! vF; 6 f 1) V COL ! 50) ,FC63 9 COL ( 63) 9 FS 5, 25 9 COL J 79) 9 F C 5 9 3 5 9 COL I 90 5 9 
F; 5,3)9C0L! iOI )9F«‘>9i5 «COLi 1 1 3 ) ,F t ^ , 1 ) ) ; 

Q_JP = ELOiJ: 

B 00 _uo = L_0.-,'N / <0„02 * R5VTENP -a- 0o6); 

DEF_U“ = OnSATR - OG_CWN; 

1<?.20_UR = K I2C_HIX5 
GO I 0 A 0 0 ^ 

/ <: i;c '‘A V* ❖ V v: i;: j;: <c C: ❖ Ci ❖ * * 0 ❖ 0 * C: v / 

/[nIOUSVkV D I SC H A G ! i\G l^ASTE INTO RIVER / 

/ )C; * V 2j: V* »:< 5'<s v;x :;c v iC: ^;< O -;> j;: ::: I.i * V sC* j:: v« ^ ❖❖ O <•* ^ / 

/ zaU V i> <* 4: 15: i‘,i * 4: C: i;« 1/ ::: -'A ijA 

FLSE :f typf=''IO' oc^ 

y VJD^OI S: GET E G ITI U? STN v G .;NSTM , 0 1 S_Z v? R i VTE NP 9 W I OTH t, OEP VH 5 (I COL <l 2 5 5 9 P H 60 2D 
., COL< 35) 9 F{ 69 2 ; 9 COL { 45) ,F{ 692 ) 9C0L{ 57) ,F{49 1) 9 COL (65) ,F (69 1) 9 
COL ( 75) ,F ( 69 1) ) ; 

r.Fr FOiT( 3005900, ;a 20_109 1 NTENP, ZILCH) (COL (25 ) ,F(69 1) 9CCL!37S 9 F (4, 1 5 9 
COL (46) ,F( 5,3 ) 9 GOi. ( 5 7 ) , F { 4 , 1 ) t X ( 19 ) 9 A U S ) 5 
!F. N=l then D0_UP=00; 

ELSE 90_UP=D0SATR-0EF_UP; 

OOSATW = 14.65 - 0.4i*IN7£NP 0 „ 008*1 NTE MP=^’^2 - 0 . 00008* I NT£MP**3 5 '' 

L 3NGTK=; UP STS-DVJNST;--;) *5 2 80.0; 

CSA = 0.75 * 'LIOTH =.- DEPTH; 

VOLUME = CSA * length; 

F L ow= Q_u p •> ; 0 1 s_ ; * i . 545 : ; 

VFL=FLni-;/CSA ; 

rr me=( length. -VEL ) /36400.0 ; 

<120_E=K :20_10*2.303; 

•<i 2 0 _or x=( ; K i 20 _up*o_up; v ;xi 20 _e *(0 is_i * 1 . 545 ) ) ) / flow; 

<220_1D=; 12o50*yFL**0.55 /;DEDTH**1.5) ; 

K22O_e=K22C_10*2.303; 

<L = <i20_NI X*1.047**’:R) VTENP- 2 O 0 O) ; 
i< 2 = X 2 2 0 _ E * 1 o 0 : 5 5 * * ; r< : y V E : • . P - 2 0 o 0 ) ; 

<3 = K320_E * :o047**;R:y7EHP-20.0) ; 

LA_iN---BQ05/( l.C-1 .0/ : 2. 71 **; 5*!<1 20_E 5)5; 

LA20 = ( ( 3OO_UP*0_UP) * CLA_IN * DIS_I5) / FLOW; 

LA = LA20 * (0.02 * RIVTE'-'.P * 0.60); 
nFF_IN = OQS.,Vl'-DO; 

00SA7R = 14.65 - 0 . 41 *R VTEMP ■> 0 . 008*R I V TE WP**2 - 0. 00008* R I VTEViP**3 ; 
DEF = { { ; DGSATR-D0_UPS*C'_UP5-!-{DEF_£f\:*0IS_l ) J/FLOW; 

IF FORK(T) = 'RES’ THEN DO; 

U = FLOV.' * ( 86400. o/vclu.’;e) ; 

Z = 0_UP * ( 86400. 0/VOLUHE ) ; 

A = Kl * K3 * '.;; 



B = i<2 + W; 

BOO_.G: L_Di-.'V = L A_! WT T=:= ; 1 ,0/ ( 2 , 71**11 A*T IME H ) <- UZ*LAD/AJ * 

c io*-i„n/i 2 . 7i**i A*‘:5r-'c) ) 3 ; 

C = DEF_IMIT - <Z*DcF5/a; 
n = ;<: / ( A-R ) ; 

£ = LA_INir - (Z*£A)/E; 

F = io0/( 2o7i**.:s*'.‘if-;£) ; 

G = Kl/i R-a3 ; 

s = la_tn;t - (z*La5/a; 

I = I „0/i 2o71**( M*riMES 3 ; 

J = ZZB; 

K = OEF ^ ?!U*'.A)/A; 

0£F_I0: DEF_DW'''! = (C«-ID*E)3*F <■ {G*H3*I ^ IJ*K3; 

FNO ; 

ELSE OQ; 

i_oi:n = (LA - * {l„o/(2.7i**(:u*time 5 )j «■ x/ki; 

FXP'l = ioO / C2<.71**il^l*Ti;-;E5 ! ; 

EXP2 = ioO / ( 2<.7:**JK2*riM£i 3 ; 

IF K2 = K1 TKEN 00; 

L = K1 * TISE; 
r5 ~ 

P = 

DEF_D'..'M = ^ DEF -c- ^ EXPl END; 

ELSE 00; 

n£F_ow\' = (:ki*la - x3 / ;:< 2 -!<D) * ;expi - exp2) <i- nx<-.v5/K23 * 

( 1 «0-l,0/( 2« 71**1 K2*TIK£ 3 3 5 * OEF* ( 1 =0/ H 2. 7 1** U<2*T IME 3 3 C 5 END; 
END; 

IF LENGTH < (0«1*5280«0) VFIEN 00_OWN = DO_UP; 

ELSE m_DWN = CCSA7P - D£r_OWN; 

PRINT_!0:OUT SXIP FD i T ! RE aGH , TYPE , UP STH, OWNS TM , FLOW ,CSA , T I KE <, X I , K2, 

00_GP , Dn_0l-;N3 ;COL 52)^»A;3 5»C0L(939A133 »C0L { 17 3 , F (6 »2 ) »COL (2839F(6j23» 
C3L(37)5F{6,i;9CnL(58}jF{6) , COL (66 3 7Ft5(,23»C0L(79)9F115j38 » COL 190 3 , 

F< 5, 3 ) 9 CQL( 101 3 vF{ ^9 1 ) jCOLI 113 3 .r 13 3 ; 
q_UP=FLOW; 

RD0_UP = L_nWN / e0o02 * RIVTENP <- 0o6); 

OEF_UP = DOSATR - D0_0UN; 

<I20_lJP = K120_NiX; 

GO TO 400". 

END; 

/ j;; i:j * V C: <c V Vf V* Ai V' •'t: i:-' * V V* '•;£ i;: 1 : 1 Jj: 2C« j:: ^ * <5 * ❖ ❖ ❖ ☆ ❖ ❖ ❖ 5^ ❖ C: CJ * O / 

/ V >;= V I-/ ^ C: V '4^ ili O i;s V J*.* ■"C'' V <: i:? 1 :: O « Cf Ci Cf <» * <1 / 

J 1 ,; vx V V V 7,'^ ^ i,s 1.x 1,; INDUSTRY UI THDRAUING WATER FROM RIVER 

/•,•£«: .S4( * -i* <<***************************************************** *****/ 

/<£<£ .>>>>:!•;< vi:-^ •.;<**************••.:************************************* *******/ 

ELSE IF TV?F=’iW= THEN DOS 

IND_UTH0R.;L;GET cC:7(U?Sr:-;9DWNSTM,WTKDR9RIVTEF3P9TJIOTH,OEPTH3 {COLI25 5 9 
F( 692 3 9 Cni; 35 3 9 p( 6925 9 C 0 LI <^5 3 ,F«6,23 , COL (57 8 9 F(<ml 3 jCOLI 65 ) , F( 69 1) 9 
C0U 75 3 9 FI69 13 >S 
IF N =1 TKEN nO_U?=GQ; 

ELSE 03_UP=D0SA7R-DEF_UP; 

LENGTK=; UPSTM-OWNST:-!5*52GOoO; 

CSA = 0„75 * UIDTK * DEPTH; 

V0LU'«:E = CSA * LENGTH; 

FLnw=0_UP-H WrHDR*lo 5458 ; 

VEL=FLnw/CSA; 
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1 I ^ t_ U MG . K / V ^ ^ /C 0^>C 0 o o V 

k 12 o_‘. 1 ;:=!Ci2c_u2 „ 

;v220_ -0= ; 12..^>0 'Vcl'"-0oG? /: .5> S 

IC220_E = ;<220_10 * 2oGCSC 

1 = :< L 2 0 _ i’ ;i '• : . o : v v 2 x ? - 2 o , o ; ; 

•c 2 = : c 2 2 0 _ E S' I o 0 1 S ^ V V ' E ; V " 2 - 2 0 o 0 J ; 

K3 = ;<320_E * .o0^-7=:=s=c:;:vv2;;2-20o0i; 

L4.20 = ECO_UP ..-2Cv;/c_gp: ; 

LA = LA20 * (O 0 O 2 =- n’JVTEM? c- 0<.60J; 

onSATR = ;.Ao65 - OoS-S-R:vT£t-'.P fl- 0«003*Rl VTEMPP*2 - 0 .OOOOSS-RI VTew)P<“»3 5 
0EF = DOS;.TR-OO_U?" 

I F F jR -: ; 7 ■ = - Ac L ’ VK2M ZZ% 

v; = 2LG'.-! * i 86*'.C0«0/V'CLL..2 i ; 

2 = Q_UP S' ■; 36-'.-.00«0/VOLL;i-.c5 ; 

A = K1 ■:• X3 <- 


8Do_:v:: :/C2^ 7 i^^u<^tiked n ^ chz^lai>/a& 

C i o 0~ I o 0 2 u 7 ^ i it 9 

C = CcF^I'X'IT - vL.-02r:/3^ 

~ t\.j. 7 iA"“3/v 

E = LA^I^JIT - IZ^LC^l^Z: 

F — i.o0/’, ^o7^^*''”vLr''^lA <lE J 'i 9 

G = x:/: b-a; ; 

H = - :Z-LA;/AS 

I = loO/i 2o7I^*-C A'-Tir^Ei; 5 S 
J = Z/3^ 

K = Gt? V 

OFF_IU: 0£F_DWM (C<*CO-En*F (J^'KDS 

EMO 9 

ELSE DOS 

L_ov:v = CLA ~ : :.o/C2o7i'i-’<*(Ki*nr;Em ^ x/xis 

EXPl = i«0 / C 2o7I--C ICl-L'TirsEn ; 

FXP2 = ioO / C 2o7I^*^*i; X2^VI.XE) ) ; 

IE ;<2 = XL TKc^3 GO; 

L = XI V tixe; 

0 = :LA-xi/xi; 


■: x^v: /x:; 

^.D’JN - DEE 

E DO; 


p = 
DEE 
ELS 
DEE 


ETsJD 
I F 
ELS 

PRINV^r.i'^PUV SXI 


EXPi - 


EHDl 


o'.:iv; = ; 

X1^:^Lf. -- X} / 

CX2-;UD ^ 

C*‘ 

CeXPl - EKP25 ^ JOIOV5/K23 R 

: o0-‘. oO/ 

u 2 o / iC 2^* * 4. 

XE ns ^ 

DEE 

S'{1<.0/H2<.71*»SK2S'7IKED 33 ? EMC; 

EMGTK < 

\ 0 J ‘ ^ L C C o G if 

TXEiv: D0_ 

DWW 

= oo_UP; 

03_cv.';c 

- i,^C3A»K ^ Cc; 

E_C’.:^:; 




03_LFs,00_0WM; JC01:2: »Ai3D oCCLC95 ,A(35 sCOLJ 1L7B ,F{6„2 J.COL528}9plI6o2U 
COL ( 3 7: ,F{ 6y 2 5 yCC_C5&: ,FC6) y CGL( 68 5 y F 5 5« 2 5 ,C0L«79) »F {5 , 3 5 » COL « 90 5 » 

F; 5, 35 y CGL; :0L 5 yF{<-y 1 5 yCCL( lI3)yF24,U}; 


Q_i;p=FLDVJ; 

B00_UP = L_D'...; / C0„02 S' RtVTEf'.P -o- 0«65; 
r>EF_UP = DOSAVR - D0_0'..'i\S 
iC 5. 2 0_0P — 1 20_s, i A y 
GO TO ADO; 


E K' D 


' V V V V w ^ 


<c * :;i ❖ 4: j:! ^2 « 3C: *<£ 5J* O ❖ j;s C: Cl O Cl C C: <2 O J> / 



/ V V' V V ^ '»s V ^ •■' *.* • . * '.* *(^ *(* V* *1^ ^ V ^ ^ V V O O ^ ^ O V V ^ ^ S ' ^ ^ f 

/**,u*».-.-*^<.-<<*^tv-.;.>:; rjgoUSTRV HI THOR 4W TNG G'ATER FROM RIVER j;;!::**-.-:* *«.•*<■■»>:: 00/ 
/** REACH IN '.■.•KICH AN INDUSTRY'S U I THOR A'.JAL AND DISCHARGE POINTS ■>*<=*/ 
are close ENOUGH TO BE CONSIDERED AS ONE POINT <=*<:OC=c:*c=>s=>C!/ 

^ ^ >]r il< W jJ: X* V ^ X/“ V V V V V V V* *1* ^ S|i ijj jjs ^ >}c i,j ^ ^ ^ Xr ^ M ^ ^ C* 0 C* ^ 

/ S}{ 5^r * Jjf s;? -if. <? i:< i;* )>■ >:« :.: s.'s * :;« * xji 

ELSE IF TVPE==:WC^ THEM DO; 

INn_svTHDRV'_cnN; GET ED I VI UP STK, DVJNSTM » WTHOR , ZI LCH H COLI 25 Ufi&tZU 
COL ( 3 5 5 5 F : 6 V 2 5 p COL i <-5 3 p r { 6 , 2 ) , X « 29 ) , A { i ) 5 ; 

IF N=l THEN Dn_lJ? = DO; 

FUSE DO_UP=DOSATR-D£F_UP; 

FLO '..'= Q _U P V ( '.: V H D R * 1 o 54 5 ) ; 

DO_DV,'N = m_UP ; 

LA20 = 33n_UP * {FLOl.'/Q_U/I ; 

PRINT_IiJC: PUT SRI? E DI T I REACH , TYPE pUPSTM , DWNSTMp FLOW p DO_U? » DO_DWN J 
;COL ( ? ) p A { 33 ,001193 pA’ 33 pCOLU 7! pFJ6p 25 pC0L(28) , F I 6, 2 ) p COL I 37i) , 
F;6pi;pC3L5I013,F(4pl)pCQL(li3),FKpl}!; 

0_JP=FL0W; 

BDD_UP = LA20; 

OEF_JP=DEF_OWN; 

<i20_UP=K120_SIx; 

GO TO ADO; 

LASTjEND START; 



APPENDIX II 


DATA FOR EACH REACH 
FOR 


JUNE 1972— SEPTEMBER 1972 
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TABLE XI 

CLASSIFICATION, RIVER MILE LOCATION, AND 
CHANNEL DIMENSIONS OF WARRIOR RIVER REACHES 


Reach 

No. 

Classi- 

fication 

Upstream 
River Mile 

Downstream 
River Mile 

Average 

Width 

(ft.) 

Average 

Depth 

(ft.) 

lA 

SD 

385.00 

382.00 

1125 

54.7 

2A 

SD 

382.00 

367.50 

1270 

71.7 

3A 

SDH 

367.50 

356.00 

1250 

83.7 

4A 

R 

366.00 

365.00 

1440 

85.7 

5A 

N 

365.40 

361.50 

1140 

32.5 

6A 

SDH 

351.50 

347.50 

■ 1240 

60.5 

7A 

R 

347.50 

346.90 

1500 

76.5 

8A 

N 

346.90 

346.30 

970 

16.9 

SA 

SD 

346.30 

345. 70 

500 

22.9 

lOA 

IWC 

345.70 

345.70 

— 

“ - 

lOB 

ID 

345.70 

345.20 

520 

19.9 

iOC 

SDH 

345.20 

344.40 

520 . 

19.9 

llA 

IWC 

344.40 

344.40 

— 

- - 

IIB 

ID 

344.40 

344.39 

470 

20.4 

ilC 

ID 

344.39 

343.70 

470 

20.4 

12A 

SDH 

343.70 

343.50 

500 

19.9 

13A 

ID 

343.50 

343.00 

530 

19.4 

14A 

IW 

343.00 

342.50 

640 

18.9 

15A 

ID 

342.50 

341.60 

500 

19.9 

16A 

SDH 

341.60 

338.41 

590 

34.9 

17A 

ID 

338.41 

338.39 

470 

49.9 

18A 

R 

338.39 

338.10 

750 

51.4 

19A 

N 

338.10 

337.00 

510 

25.0 

20A 

IW 

337.00 

336.70 

410 

25.3 

2 LA 

ID 

336.70 

335.00 

350 

25.5 




TABLE XII 


INPUT DATA FOR WARRIOR RIVER REACHES 


Reach 

No. 

Mon th 

Input 

Flow 

(cfa) 

Rivei* 

Tempera- 

ture 

(°C) 

BOD 5 

Input 

(ppm) 

D. 0. 
Input 
(ppm) 

20 °C 
Input 
(base 10 ) 
(days*”^) 

Teripera-' 

ture 

Input 

(OC) 

BOD 5 

of 

River 

(ppm) 

Do Oo 

of 

River 

(ppm) 

ki 20 °C 
of River 
(base 10 ) 
(days"”^) 

lA 

June 

1895.0 

27.0 

2.5 

6.0 

0.060 

27.0 





July 

2150.0 

27.9 

3.2 

6 . 7 

0.060 

27.9 





August 

1615.0 

30.0 

1.0 

6.6 

0.060 

30.0 





September 

1975.0 

28.0 

3.0 

5.9 

0,060 

28.0 




2A 

June 

65.0 

27.7 

2,5 

7.0 

0.060 

25.6 





July 

100.0 

29,0 

3.5 

5.2 

0,060 

28.9 





Augus t 

100.0 

30.0 

3,5 

7.1 

0.060 

28.9 





September 

125.0 

28.5 

2.5 

6.5 

0.060 

25.0 




3A 

June 

10.0 

28.5 


8.0 


20.0 





July 

10.0 

28.0 


8.0 


22.2 





August 

15.0 

29.1 


8.0 


22.8 





September 

10.0 

28.9 


8.0 


22.8 




lik 

June 

■ 1970.0 

28.5 





1.0 

3.7 

0.060 


July ‘ 

2260.0 

27.5 





1.5 

3.4 

0.060 


August 

1730.0 

29.1 





1.5 

5.1 

0,060 


September 

2110.0 

28.9 





1.5 

3.0 • 

0.060 


00 

K> 



TABLE XII (Continued) 


Reacli 

No. 

Month 

Input 

Flow 

(cfs) 

River 

Tempera- 

ture 

(°C) 

BOD 5 

Input 

(ppm) 

D. 0. 
Input 
(ppm) 

5A 

June 

0.0 

27.3 




July 

0.0 

26.2 




August 

0.0 

28.2 




September 

0.0 

27.8 



6 A 

June 

15.0 

26.7 


5.0 


July 

25.0 

26.0 


5.0 


August 

10.0 

29.0 


4.6 


September 

15.0 

28.3 


4.8 

Ik 

June 

1880.0 

28.3 




July 

2270.0 

27.5 




August 

1730.0 

28.9 




September 

2040.0 

28.0, 



8 A 

June 

0.0 

25.7 




July 

0.0 

26.7 




August 

0.0 

27.2 




September 

0.0 

27.2 



9A 

June 

25.0 

25.7 

2.0 

8.0 


July 

35.0 

26.7 

2.5 

8.0 


August 

25.0 

27.2 

2.5 

8.0 


September 

15.0 

27.2 

2.6 

9.4 


k3_ 20 °C 

Tempera- 

BOD 5 

D, 0. 

20 

Input 

ture 

Of 

of 

of River 

(base 10 ) 

Input 

River 

River 

(base 10 ) 

(days"“^) 

(oc) 

(ppm) 

(ppm) 

(dc.ys"-^) 


0.055 

0.055 

0.055 

0.055 

22.2 

23.3 

26.7 

26.7 


2.7 

4.6 

0.095 

3.5 

4.3 

0.095 

4.0 

3.9 

0.095 

4.0 

4.4 

0.095 


0.065 

0.065 

0.065 

0.065 


0.060 

25.6 

0.060 

28.9 

0.060 

28.9 

0.060 

25.6 


CO 

u> 



TABLE XII (Continued) 





River 





Input 

Tempera- 

BOD 5 

D, 0. 

Reach 


Flow 

ture 

Input 

Input 

No, 

Month 

(cfs) 

(°C) 

(ppm) 

(ppm) 

lOA 

June 

- 0.43 





July 

- 0.43 





August 

- 0.43 





September 

- 0.43 




lOB 

June 

0.40 

26.5 

610.0 

0.0 


July 

0.43 

27.3 

450.0 

0.0 


August 

0.43 

27.8 

280.0 

0.0 


September 

0.42 

27.2 

410.0 

0.0 

IOC 

June 

0.20 

26.5 

1550.0 

0.0 


July 

0.20 

27.3 

1550.0 

0.0 


August 

0.20 

27.8 

1550.0 

0.0 


September 

0.20 

27.2 

1550.0 

0.0 

llA 

June 

- 1.12 





July 

- 1.12 





August 

- 1.12 





September 

- 1.12 




IIB 

June 

0.21 

26.5 

5.0 

3.0 


July 

0.21 

27.3 

5.0 

3.0 


August 

0.21 

27.8 

5.0 

3.0 


September 

0.21 

27.2 

5.0 

3.0 


ki 20 °C 

Tempera- 

BOD 5 

Dc 0. 

ki 20 °C 

Input 

ture 

of 

of 

of River 

(base 10 ) 

Input 

River 

River 

(base 10 ) 

(days“l) 

(°C) 

(ppm) 

(ppm) 

(days*~l) 


0.005 

33.0 

0.005 

33.0 

0.005 

33.0 

0.005 

33.0 

0.025 

33.0 

0.025 

33.0 

0.025 

33.0 

0.025 

33.0 


0 . 

.025 

34.4 

0 , 

,025 

34.4 

0 , 

.025 

34.4 

0 , 

,025 

34.4 



TABLE XII (Continued) 



River 



ki 20 °C 

Tempera- 

BODr 

D. 0. 

20 °C 

Input 

Tempera-’ 

BOD 5 

D. 0. 

Input 

ture 

of 

of 

of River 

Reach Flow 

ture 

Input 

Input 

(base 10 ) 

Input 

River 

River 

(base 10 ) 

No. Month (cfs) 

(°C) 

(ppm) 

(ppm) 

(days“^) 

(°C) 

(ppm) 

(ppm) 

( ' ys"^) 


lie 

June 

0 . 70 

26.5 

84.0 

0.0 

0.085 

33.0 


July 

0.80 

27.3 

52.0 

0.0 

0.085 

33.0 


Augus 1 

0.78 

27.8 

54.0 

0.0 

0.085 

34.0 


September 

0.72 

27.2 

46.0 

0.4 

0.085 

33.0 

12A 

June 

70.0 

26.5 


7,5 


25.0 


July 

110.0 

28.0 


7.0 


28.9 


August 

110.0 

27.8 


7.0 


27.8 


September 

140.0 

27.2 


7,5 


24.4 

13A 

June 

3.50 

26.5 

2.2 

5.0 

0.002 

24.0 


July 

3.50 

28.0 

2.2 

5.0 

0,002 

24.0 


August 

3.50 

27.8 

2.2 

5.0 

0.002 

24.0 


September 

3.50 

27.2 

2.2 

5.0 

0.002 

24.0 

14A 

June 

- 18.00 

26.5 






July 

- 18.00 

28.0 






August 

- 18.00 

27.8 






September 

- 18.00 

27.2 





15A 

June 

12,00 

27.0 

100.0 

0.0 

0.135 ’ 

35.0 


July 

12.00 

27.5 

100.0 

0.0 

0.135 

35.0 


August 

12.00 

27.8 

100.0 

0.0 

0.135 

35.0 


September 

12.00 

27.2 

100.0 

0.0 

0.135 

35.0 



TABLE XII (Continued) 


Reach 

No, 

Mon tli 

Input 

Flow 

(cfs) 

River 

Tempera- 

ture. 

(°C) 

BOD 5 

Input 

(ppm) 

Do Oo 

Input 

(ppm) 

20 °C 
Input 
(base 10 ) 
(days~^) 

Tempera- 

ture 

Input 

(°C) 

BOD 5 

of 

River 

(ppm) 

Dc 0 . 

of 

River 

(ppm) 

20 °C 
of River 
(base 10 ) 
(days”^) 

16A 

June 

0.6 

27.0 


4,0 


21.0 





July 

0.4 

27.5 


4.0 


23.0 





August 

0.5 

27,8 


4.0 


22.0 





September 

0.6 

27.2 


4.0 


21.0 




17A 

June 

0.32 

26.7 

58.0 

7.0 

0.100 

27.0 





July 

0.50 

28.5 

40.0 

8.0 

0.100 

28.0 





August 

0.38 

27.8 

28,0 

9.0 

0.100 

28.0 





September 

0.50 

27.2 

60.0 

7.0 

0.100 

27.0 




18A 

June 

2150.0 

26.7 





2.5 

6.6 

0.040 


July 

2630.0 

28.5 





1.5 

5.7 

0.040 


August 

2000.0 

27.8 





1.0 

5.4 

0.040 


September 

2320.0 

27.2 





2.0 

6.7 

0.040 

19 A 

June 

0.0 

26,7 







0.020 


July 

0.0 

28.5 







0.020 


August 

0.0 

27,8 







0.020 


September 

0.0 

27.8 







0.020 

' 20A 

June 

- 0.40 

26.7 









July 

- 0.40 

28.5 









August 

- 0.40 

27.8 









September 

- 0.40 

27.8 










TABLE XII (Continued) 


Reach 
No . 

Month 

Input 

Flow 

(cfs) 

River 
Tc ipera- 
ture 
(°C) 

BOD 5 

Input 

(ppm) 

D. 0. 
Input 
(ppm) 

21A 

June 

0.46 

26.7 

22.0 

2.9 


July 

0.46 

28.5 

22.0 

2,9 


August 

0.45 

27.8 

27.0 

3.0 


September 

0.46 

27.8 

27.0 

3.0 


20 °C 

Tempera- 

BOD 5 

D. 0 . 

kj^ 20 °C 

Input 

ture 

of 

of 

of River 

(base 10 ) 

Inout 

River 

River 

(base 10 ) 

(days”^) 

(°C) 

(ppm) 

(ppm) 

(days'll) 


0.025 35.0 


0.025 

35.0 

0.025 

35.0 

0.025 

35.0 


00 



APPENDIX III 


DETAILED RAPS OF THE SECTION OF 
THE WARRIOR RIVER USED IN THE 


MODELING STUDY 
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